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4, study is made of compressible perfect fluid motion in 
turbomachines having infinitely many blades and a general theory 
is developed. Ап underlying concept of the theory is that force 
fields which represent the action of infinitely many blades belong 
to a special ciass described as pseudo-conservative and can be 
expressed as the product of a scalar function and the gradient of 
a potential. ‘he scalar function is simply the rate at which energy 
is imported to the fluid by the blades, and the potential is simply 
the family of the equations for the blade surfaces. The introcuc- 
tion of these two functions to express the force field casts an en- 
tirely new light on problerns of mixed-flow turbomachines having 
infinitely many blades of arbitrary shape. 

In the formulation of the problem the non-linear action of 
rotationality and compressibility is regarded as a force tending 
to displace the streamsurfaces from their irrotational, incom- 
pressible position. It is shown that the character of the problem 
is determined by a governing velocity: the velocity relative to the 
blades where blades are present, or the meridional velocity, 
where blades are not presen:. Where the governing velocity is 
subsonic the problem is essentially ellíptic, where supersonic, 


hyperbolic. 
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The theory and the examples lead to conclusions which 
ат Й 


are believed to explain in part the unexpected efficiencies ob- 


served for compressors having transonic governing velocities. 
These conclusions, which indicate that transonic compressors 
could perhaps be profitably developed, are as follows: The de- 
flection of the streamsurfaces induced by a given strength of 
vorticity at a certain point in the flow has one sense when the 
governing velocity at the point is subsonic, the opposite sense 
when it is supersonic, and becomes zero as it becomes sonic. 
The deflection of the streamsurfaces brought about by a given 
distribution of vorticity in a region is less when the governing 
veiocity in the region is transonic than when it is entirely sub- 
sonic or entirely supersonic. 

Юхагар!е5 of incompressible flow through a mixed flow 
compressor with prescribed blades, and subsonic and transonic 
flow through actuator disks, were solved by the method of 
finite differences, applying simultaneously the relaxation 


technique and an iteration process. 
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1. INTRODUC MEN 


In a real turbomachine a viscous, compressible fluid fiows 
through an axially symmetricai channel. Inc region of this channe! 

a systern of blades, either stationary or rotating about the uxis of 
symmetry, acts on the fluid. The field of the ilow, being bounded 
py the surfaces of the blades as well as by the channel boundar.es, 
is not circumferentially uniform. The total energy of the fluid may 
vary írom point to point in the field, the flow is genera-ly rctational, 
and if the blade system rotates, the flow is unsteady. Whe real 
problem is thus a very formidable one, and cannot be approachec 
exactly by methods known today. tf viscous effects are neglected 
the problem is greatly simplified. because of this, and because 
viscous effects actually are very slight except near the boundaries, 
the assumption of zero viscosity has always been made. 

/. brief but excellent review of the earlier important investi- 
gations of the flow in turbomachines is given by Marble (2eference 1). 
Zarlier investigations were concerned primarily with the flow through 
à typical annulus of small radial extent and hence treated the flow as 
essentially two-dimensional, e. р. Traupel (Reference 2). Inierfer- 
ence between the flow in neighboring annular regions was assumed 
negligible, a condition which is fulfilled only if the centrifugal ior 


and the radial pressure gradient forces are in balance, tnat is, if 
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the flow is that of a vortex centered on the axis of rotation. This 





is a rather severe restricticn, too strong to be irnpased in most ap- 
plications, - 

The first detailed analysis of the three-dimensional incom- 
pressible flow in turbornachines was given by Neyer (Reference 3). 
An exact solution was obtained for machines in which the flow is ir- 
rotationalupstream and downstream of the rotor, although the flow 
within the rotor rnay be rotational. h eyer achieved circumferentia! 
uniformity by assuming an infinite number of very thin blades, and 
then modified this result by a Fourier analysis to obtain the solution 
for a finite aumber of blades. This wes an important contribution 
but the extension to arbitrary blade shape, with rotational flow down- 
stream of the rotor, has not been made. , 

Marble (References 1 and 4) linearized the rotationai incom- 
pressible problern with an infinite number of blades by the assump- 
tion that the vorticity is transported aiong the streamlines of an ir- 
rotational flow within the same boundaries. He provided examples 
of axial flow and conical flow, and stated that the simple linearized 
solution was sufficiently accurate if the vorticity effects were not 
large. A second order linearization was given to handle problems 
in which the vorticity effects are large. Using the simple linear- 


izations Marble considered three interesting problems: the mutual 


interference of neighboring blade rows ina multistage axial turbo- 
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'bitrary. Tbe two assumptions Which snake 
the analysis possible are: the fluid considered is non-viscous, and 
the number of blades is infinite. The viscosity of the fluid is unim- 
portant except within the boundary layers along all the bounding sur- 
faces, and these boundary layers are thin if pressure gradients are 
favorable. furthermore, an elementary consideration of boundary 
layers often presupposes a knowledge of the velocity distributions 
such as obtained by this analysis. The assumption that the number 
of blades is infinite makes the flow field circumferentially uniform. 
This uniformity makes the solution easier in that the flow depends 
only on two coordinates, oui at the same timo it introduces arti- 
ficiai complications. Generally there wiil be a discontinuity in 
the flow as it enters this blade region if the number of blades is 
infinite. Furthermore the force o! the blades on the fluid is not 
applied on individual blade surfaces, but is distributed throughout 
the blade region and actually acts as a body force field distributed 
throughout the region. The circumierentially uniform flow field 
may oe considered as the limiting case of a flow acted upon by a 
very large number of closely spaced blades of negligible thickness. 

Two distinct problems occur in practise. "he design prob- 
lem arises wher it is desired to design a machine for a particular 


purpose. lhe analysis problem arises when it is desired to 














investigate a given machine under given operating conditions. These 
problems were classified by Marble (Reference 4) as the indirect and 
the direct problem, respectively, in analogy with the three-dimen- 
sional wing theory. This classification is convenient and very sig- 
nificant. Here any problem in which the blade shape is prescribed 
will be classified as a direct problem, any other problem as an in- 
direct problem. The full significance of this distinction becomes 
clear when the circumstances under which the idealized problem is 
actually cornparable to a real problem are understood. 

In the real machine, which has a finite number of blades, 
each blade transmits a force to the fluid by maintaining a discontin- 
uity in pressure across iis two surfaces. If no viscous force is 
present this blade force must be normal to the blade surface. In 
idealizing the problem these concentrated blade forces are essen- 
tially replaced by a body force field. Clearly this body force field 
must be normal to the family of blade surfaces throughout the region 
of the blades if the two problems are to be comparable. It is im- 
portant to note that this imposes a purely geomctrical restriction 
on the body force field. It will be shown that the necessary and suf- 
ficient condition that it be possible to construct a family of surfaces 
which are everywhere perpendicular to a given force field is that the 
force field is everywhere perpendicular to its own curl. lí the blade 


shape is prescribed, as in the direct problem, this offers no difficulty. 











the things tc be determined, the prescribed quantities, which for 
example, may be components ot the force field or the energy dis- 
tribution, must be prescribed in such a way that the existence of 
the family of blade surfaces is assured beforehand. This necessity* 
was noted by lLieyer (Reference 3) in connection with the srescrip- 
tion of vorticity in the indirect problem, but was overiccked by 
Reissner (Reference 6), Gravalos {Reference 7), and several others. 
When the streamwise extent of the working region is m.uch less than 
its transverse extent, that is, if the Siade aspect ratio is large, as 
іп rnost exial flow problems, the idealized flow will be cornvarable 
to the real flow even thouga the blade existence has not been assured. 
Thus Marble (Reference 4) was able to obtain good approximate 
solutions of the indirect axial flow problem: by prescribing blade 
force components. However in mixed-flow probieras the aspect 
ratio is not large, and in the indirect mixed-flow problem the exist- 
ence of the blades must be assured beforehand. 

In the rea! problem one boundary condition is that the fluid 
does not flow through the blade surfaces. This cannot be applied 
as a boundary condition in the idealized problem since the blade 


surfaces, being infinite in number, are not boundaries of the flow. 


*Informaticn lately received reveals that this condition for integra- 
bility of the blade surfaces was stated by Pauersfeld as early as 
1905 in a paper (Reference 10) not available at tne time of writing. 
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ment on the Glade surfaces. This constraint, quite different from 
the boundary condition of the rea? problem, plays an important part 
n the mathematical formulation of the ІМегііте protles.. The 
geoiretricai relation which expresses the contiraist ig equal ir 
significance iu the cther equations end must Бе соте габ візі. 
taneously with thems. In addition, with uniform inlet conditions and 
an incompressible fluid, the differential equation for the flow 15 
Nagar in the region of the blades, alcnough it is genera iv non-linsar 
downstre in of this region. This fact sas tied? eyer (Referenca 3) 
tc Над ехасі solutions for lows iv which Che vorticity Fawnstncam 
of the biades ia zero although the flow сау Se rotational and vary 
comp icated ia the region of the bises. 


Го suci5rigse, the ideajlicad problera must satisfy the € QUE - 


equetici, the wentropic prosauco-ensit ty Termtion, ang theese boun- 
diry conditions cf the rea! problem: not applied on the blade surfaces. 
The two additional requirements which must be fulfilled if the $2. 

lution of the idealized problern is to be co rpe?rzble te the real prob- 
6сп аге: (ће wody force field ous! be ne emal to the family ئم‎ Slade 
Surfaces; and the velocity vector field срив: be ocrrallel to the farm. 


ily GY Glade surfaces. 
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supplementary conditions can easily be imposed. Sor the indirect 
problem in which force components are prescribed these require- 
ments should be expressed in terrns of the force components as 
follows: the body force field must be normal to the velocity vector 
field; and the force vector must be normal to its own curl. These 
assure the existence of a family of surfaces which can be chosen 
to fulfill the two requirements stated above for the direct problem. 
If in the incirect problem cther quantities such as energy distribu- 
tion are prescribed, equivalent requirements must be satisfied. 

The idealized problem is formulated in terms of a stream 
function for the velocities in the meridional plene, and the result- 
lug differentiai equation for the stream: function is written as a 
second order non-homogeneous partial differential equation by re- 
garding the non-linezr terms, representing the effects of rotationality 
and compressibility, as forcing functions tending to displace the 
streamsuriaces. The differential equation is then replaced by a 
finite difference equation which is solved by a simultaneous appli- 
cation of the relaxation technique cf Southwell (Reference 8) and an 
iteration process. 

The general theory is developed in Parts Ii through VI, and 
the difference formulation and examples are presented in Parts 
Vil, VIII, IX, and X. The exainples were conceived in order to 


demonstrate separately the different phases of the problem. Part 








„б 
VIII is concerned with incompressible flow through a region of 
prescribed blades and is primarily an exernple of rotationality 
in a mixed-flow compressor. Part ix is offered as an example 
of the effects cf compressibility when the vorticity distribution 
is approximately constant. Part X is an example of the inter- 
action of rotationality and compressibility efíects when tne total 


velocity is transonic but the governing velocity is subsonic. 
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II. CCCRDINATES, NOTATION, SYMROÜLS 












The flow is described (Figure 1) for the most part in an ab- 
solute cylindrical coordinate system r, 8, z by the velocity con.pon- 


ents u, v, w respectively. The absolute velocity vector is: 


ЖЕУ ғу”, 
where /, /g, and /z are unit vectors in the r,6, and z directions 


respectively. The vorticity components for circumferentially uni- 


form flow are respectively: 


ES | _ ди dw sow Rv (1 
ЖАЛЫ ри = F JF ) 


and the vorticity vector, A =vxV, is: 


О = rS +1 + Tz | 
4,6 





Figure 1 

The Velocity and the Vorticity Cornponents in the Absolute Coordinate 
System 

The scalar functions u, v, w, § »1,§ are independent of the Ө 


coordinate in accordance with the requirement of circùmferential uni- 


formity. 











The forca of the blades on the fluid is given by 


ITA E 


lengt! 
and has the dimensions of force/unit mass, or اس‎ 


time* 
At times the reiative cylindrical coordinates г, $, z, where 





$ is measured with respect to the rotor, will be used. If the angular 
velocity of the rotor is w , the absolute and relative velocities are 
related by : 

do -dé = dt (2) 
The velocity and vorticity vectors in the relative system are respec- 
tively: 


W = V-iywr = T, u siglv-wrl)rizw , and 


fl # {-Һго = 4 #7) # % (4-2‹0) 
The intrinsic coordinates n, 8, s will often permit a more 
concise presentation (Figure 2). nis distance normal to the merid- 
ional streamline and s is distance along the m:ridional streamline. 


The velocity and vorticity vectors are respectively: 





VE Ti 5 сањв and 
iz 205 + dE ED 8 572, , 
where 
Х2, = = pe, ay = ДЕК, 
2--6%% - E " 


and where Ke is the curvature of the meridional streamline. 


The family of the blade surfaces is given by a scalar function 








Meridional Trace 
\ of Streamsurface 





Figure 2 


The Velocity and the Vorticity Components in the Intrinsic Coordinate 
System 


oí position, called the blade surface function: 

/8(ғө,2) E (г,2) (4) 
A blade surface is defined by /3 = constant. The function Қт ,2)- (9-/8) 
defines the family of the traces of the blades іп a meridional plane. 


The normal to the blade surface is 


TB = -д o - 2% (5) 
Since the blade force F is always norrnal to the instantaneous sur- 
face of the blades it can be expressed as 

E Ar, Z] W rz) (6) 
This defines the scalar function A(r,z). From the second equation 
of motion it will be shown later that ФЛ is the local rate at which 


energy is added to the fluid by the action of the blade forces. /3 





-1 
and f are non-dimensional, 748 has dimensions of (length) ‚ апа 








E 


. | | В length* 
143 dimension 21 а. 
time“ 





Note that vS , tne normal to the blade suría ce ‚is nots unit 
vector. Instead it is chosen so that its tangential comiponeat is amit. 
Thus the magnitude of the tangential component of the force vects 
is simply Fe = HA. 


t.iscellaneous - vınbois: 


a 





M 
t 
т 


y 


y 
n 


М 


Z е mies 
2 


p = static pressure 


р = dengity 

2 = speed of sound 

' = ratio of specific heats 
t = time 


E = Lv” Je = total energy of the fluid 


( n subscript denotes condition far upetream of the region 
of the bDiades, or at the upstream boundary 


( subscript denotes condition at the trailing edge of the 
blades 


( ).( ) , partial differentiation with respect to r and 2 ro- 
збасбіугіу, unlesy otherwise defined 


а, = angle between the blade surface and the ¡meridional 
plane measured ina plane where z is constant 


a, = angle between the blade surface and the meridional 


Е 
plane ¿measured in а plane tanwent to the surface iere 
r i£ constant 
A, = augle betweeu the blade surface uni the meridional 


pi@ne Messured in a plane asrınal to their line of inter- 
section 
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curvature of the meridiemai Stream line 


ti 


H 
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curvature of the normel to the mieridianel streamline 


meridionai velocity component 


z'"residual" of the difference ooerztcr 


- differeuce operator 


circulation about an annular vortex ring 


angle between meridional velocity and axis 


influence coefficients 
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Mach number 


outer radius 


inner radius 


rr; к 
= L = non-dimensional radius 
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ІП, ГИЯ PS:UDO- CONSERVE TIVE FORCE rio 


|ND THE BLADE SURES.CZ FUNCTION 


A) The Necessary and sufíicient Condition for the Existence of the 
Family of Blade Surfaces 











з 


The body force field Y is normal tu the family of biade sur- 


мазы 


ces 2 , hence y and vS are parallel and can be related by 


к" 
р, 


Flr,2) = A(r,2) 781r,8,2) (6) 
where A(r,z) is a scalar function, and v/3 is a function of only r 


and 2 because of circumferential uniformity. writing 


and taxing the curl of both cides there results: 
(x) АҒ + я? (хх | 3 
Then forming the scaler products with > 
E а у= 
and since ДА is not zero unless |F] is zero. 
F рх F по (2 (7) 
This is the necessary condition for the existence of the family of 
blade surfaces and is the additional relation which must be fulfilled 
in the indirect probiem if the idealized problem is to be comparable 
to a real probiein. 
[n order to discuss the question of blade existence suppose 
that in the indirect problem force components are to be prescribed. 


Tirst consider flow through the channel with no force present. The 


flow is completely determined by the equations of motion, the con- 


* The applicability of Equation (6) and the subsequent derivation of 
Equation (7) were pointed out by Dr. Arthur Zrdélyi in personal 
conversation. 
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tinuity equation, a pressure-density relation, and the necessary 
boundary conditions. Now consider flow through the same channel 
with forces present. Three unknowns, the three force components, 
have been added to the vroblem. One equation, the requirement 
that the relative streamline be normal to the force field, has been 
added: 

D (8) 
at this point it appears that two components of the force can be 
prescribed arbitrarily. The third component would then be expressed 
by means of Equation (8), and the problem would be completely de- 
termined as in the channel with no force acting. The solution would 
describe a flow under the action of the prescribed forces, but would 
not necessariiy be comparable to a flow acted upon by blades because 
the existence cf the blade surfaces has not been assured. If Equation 
(7) is imposed the problem is redundant, for with two force compon- 
ents prescribed there are more equations than unknowns. An example 
will give physical meaning to this difficulty. Suppose the two com- 


ponents P, and F, are prescribed as: 


9 
F- = 4% p "n = ИЕ (Erz) 

F, can be expressed in terms of these and the velocity by Squation 

(8) as: 


2. _ VE (9). 








ТЕ 





The three force components are then known and the problem is com- 
pletely prescribed as in the channel with no force acting. The solu- 


tion will give u, v, was functions of r and z. Now the equation as to 


the existence of blade surfaces is raised. quation (7) for this case is | 





зе б = 0, 
from which 
^7. 2 = . Ғы 
m = function of z only (10) 


The problem is clearly over prescribed, for F, is restricted when 
Equation (10) is imposed. Furthermore, it is clear what should be 


с FZ з , 
із zero, then —— can be a prescribed function 


prescribed. ПГ. БЕ 
0 


of z only. From Squations (9) and (10), 


vou Tr Кт 1 
——. = - — = -r x function of z 
ум FO 


This function oí z that may be prescribed represents an angle be- 
tween the relative velocity and the rneridional plane and actua ly is 
the biade shape parameter Е (2) defined below. Thus when the exist- 
ence of blades is assured the indirect problem with force components 
properly prescribed is equivalent to the direct problem with the 
blade shape prescribed. 
B) The Most General Form of the Body Force Field 

The reason that the two force components cannot be indepen- 
dently prescribed is that the force field, although non-conservative, 
is restricted in form. The most general force field that will satisfy 


Equation (7) is given by Equation (65). This is the most general force 





«ЛЕ 
field for which a family of blade surfaces actuaily exists and is 
therefore the most general force field possible in turbomachine 
problems where the number of blades is infinite. lor this reason 
ii is provosea for ell such problems that the body force field be 
expressed as the product of a scalar function and the gradient of 
a potential: 


LIE eee Oo Ы (11) 


Tris force field hasa special piace Lbatween conservative fields 


ни 


and non-conservative fields. .. conservative force field can be 
expressed in terıns of one scalar function, its potential. A non- 
conservative force field requires three scalar functions for expres- 
sion, one for each cornponent. ¡he force fielc of ihe turbomachine, 
being intermediate be,ween these in that it requires two scalar 
functions for expression, might be called a paeudo-conservative 
iorce fieid. Actually its vectors have the same direction as a con- 
servative field, but not the same magnitude. 

This fact that the force field depends on only two scalar 
functions indicates why, in theindirect problem with force compon- 
ents prescribed, tne two components cannot be prescribed indepen- 
dently. It aiso indicates that considerable simplification might 
result from considering the functions À and f instead oi the three 
force components. 7nis is actually tne case. The introduction of 


tnese two functions for the force vector field casts an entirely полу 


light on turbomachine probiems in which the number of blades is 
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infinite. Wer example, the heretofore more difficult direct prob- 
lem becomes relatively simple. MEL = tne functions A and 
Е пауе special physical significance: DA(r,2) represents the local 
time rate at which energy is added to the fluid by the. Ration оғ го- 
tating blades, and f(r ,Z) describes the blade shape completely. 

irom Equation (11) the force field is finally expressed as и 
' 
ER VD EN | (12) 
V; hen the force field is expressed in terms of the functions 
A and íf the question of biade eridtence does not arise and the direct 
and indirect problems becorne quite ciear. / direct problem is one 
in which the blade shape f(r,z) is prescribed and A(r,z) is an un- 
known determined by the second equaticn of motion after the solution 
is complete. The function A does not appear in the equation for the 
_ 
meridional flow. An indirect probleim may be one in which the local 
rate of energy input, A (r,z) is prescribed, and the blade surface 
function, f(r,2), is an unknown dependent variable. 
C) The Fhysicai Meaning of the Blade Surface function 4 ; 
It is of interest to note the geometrical significance of ,3 
and f and their derivatives. The absolute value of 7 is 
rel s X d o» Rt. zu (13) 
where Әсе апу C -rf. 


The relation between the vector 7/3 and the blade surface is shown 


in Figure 3. Using the definitions of ар, а,, а аз given in Part 


m 
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lI it easily follows that 


QU = tan a. 


(1é) 


N 
it 


E “tan - 
do 


with the sign convention shown іп гір. 3. 






Trace of 


/ Blade 


- (-f,) 





Looking Downstream onto wooking Radially inward onto a 
a Plane z = constant Plane Tangent to a Surface 
r = constant 


Figure 3 
The Blade Surface Function f(r,z) 


In the same way 
Ж; “2 - “әр” ох, (15) 
and from Equation (13) 
г/Р|8/ = sec. *n (16) 
Thus R and Z are merely the tangents of certain angles be- 
tween the blade surface and the meridional plane. 
The functions ( and f were defined in such a way that 6 = 
constant is the equation of a blade surface and f = constant is the 


equation of the trace of a blade surface in the meridional plane. 





E. 





If the blades are so-called radial blades , that is, if the blade sur- 
faces are generated by lines that are normal to the axis ‚ the function 


f is independent oí r. in this case f, = 0 and f, is a function of z 


only. 




















iV. THE HYDROUDYNAMICAL EGQUATICNS 












-) The Equations in the Various Coordinate Systems 
The hydrodynamical equations for the steady, adiabatic, 
axially symmetrical, and circumferentially uniform flow of a non- 
viscous fluid acted upon by a pseudo-conservative* body force field 
are given in the various coordinate systems previously described, 
Vector Forms: 
The equation of motion: 

ИРИ = - жур + Ду (17г) 

/ 42 = "LUV e тр - Ар/8 (17b) 
The continuity equation: 

8 2y =o (18) 
The isentropy condition: 

(5) -0: a?s $2 - (3-1) | £ с: (19) 
Absolute Cylindricai @ordinates: 


The equations oí motion: 





pou ОС... کے‎ 20a 
быел уу = РОЗА АЛ, (20а) 
дгу 2 ги = 
Urn 5 A (20b) 
“эм F WS = P dz > (20c) 


The continuity equation: 
mot) E (21) 

Force‏ و ےک 

ЖА pseudo-conservative snes field is defined as one which can be 

expressed as the product of a scalar function and the gradient of г 

potential. See Part III. 




















Aosolute Intrinsic (low) Coordinates: 


m 


The equstions of motion: 


9; bo - sinc --£$É-A5 (22a) 

A (22b)‏ = 2 وو 

Ks 2 ЕА = -£ $5 - Af, (22e) 
The continuity equation: 

So, + Pas Ky b o espe = о 2 (23) 


B) An Integral of the Equation of Motion 


wince the flow is steady the equation of moticn can be inte- 





grated along each streamline. forming the scalar product of V and 


Equation (17b), a scalar equation is obtained: 


ИРИНУ ую Дуга = О (24) 
From Zquations (6) and (8) the condition that the relative velocity 


vector and the force vector are perpendicular is 


W: B= (7-7 or) -v 4 ="0 (25) 
Consequently 
Ив < 22, and v-wr = rfiutrhw (26) 


From the second equation of motion (Equation (20b) ): 





ше и 52 pw 2 = И Ру(Ри) (27) 


Since the flow is isentropic the density is a function of pressure only, 


р = р(р), and the pressure term of Equation (24) is: 


Шу» - х/ЯЕ (28) 
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DE e 
Using Equations (26), (27), and (28), Equation (24) becomes 
У тіж у'я | 98 - Drv] = o (29) 
Integrating along the meridional strearnlines an energy equa- 


tion analogous to the Bernoulli equation is obtained. 





+ pue. ТЕ O FY = constant on each streamsurface (30) 
Upstream and downstream of the blades rv is constant on each stream- 
surface. The total energy is E = 44%; 22 . The term z v* is the 
kinetic energy per unit mass and SE is the enthalpy. The change 
of wrv which occurs along a streamsurface in the region where hg 

lades act is the change in the total energy cí the fluid occurring on 
the same streamsuríace. 

The constant in Equation (30) will have the sarne value on all 
Strearnsurfaces if the total energy is uniform and the meridional com- 
ponent of vorticity is zero upstream of the blades. Under these con- 
ditions Equation (30) becomes 

- vês fL- дру = constant (31) 

Ina more general case the flow may be rotational and may 
have a non-uniform energy distribution at the upstream boundary, 
the inlet. Using the subscript ( ) to denote these inlet conditions 
the "constants" in Equation (30) can be evaluated. 

= Ves [EP - Фги + £, - Ә/үғи, (32) 
The functions £, and w (rv) are Constant on еасһ streamsurface and 


generally each has a different value on each streamsurface. They 











are known. 
The gradient of Equation (32) is taken in order to eliminate 
the total energy gradient, v7 E, from the equation of motion: 
vE = ТЕ 7h, - Polen), (33) 
Since the functions Z, and (rv), are constants on each streamsurface 
their gradients are normal to the streamsurface and Equation (33) 


can be written as 


= А, = / 
PE = varv + i, = - 2 ofr (34) 


where n is distance normal to the streamsurface. 
. C) Two Identities 
Two identities which can be proved for axially syr:metrícal 
flow using only the definition of vorticity are: 
(27 о; у= О (35) 


„ern (36) 


Рогу 

lhe vector j, ar has the dimension of velocity and acts 
tangentially. 
D) The Concept of Free and Bound Vorticity 

The concept of free and bound vorticity provides a useful 
means of discussing the properties of rotational flow. Although ro- 
tational flow problems are generally non-linear, they may be linear 
if the vorticity is bound. The non-linearity ті physically from 


the dependence of the solution upon the distribution of vorticity 
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which is determined, in turn, by the manner in wnich the vorticity 
is transported by the velocities of the solution. Some knowledge of 
the vorticity field can be obtained üirectly from the equations of 
motion, 

Consider first rotational flow in a region where no forces 
act on the fluid. The equation of motien is simply (from Equation 
(17b) }: 

Их = ve V2 +S vp= ve (37) 


U the total energy is uniform on the upstream boundary it is uniform 


everywhere in the field and vi > 0. This is the simplest example 
of free vorticity, one in which the vorticity vector and the velocity 
vectors are parallel*. If the total energy of the fluid is not uniform, 
the velocity and vorticity vectors are not parallel. They are, how- 
ever, perpendicular to the gradient of the total energy. 

Consider next rotational flow in a region where 4 pseudo- 
conservative force field acts on the fluid. The equation of motion 
is (from Zquation (17b) ): 

ИхО < vE-Av,8 (38) 
With Equation (33) this becomes 
V*f2-vurv-v 0 (ги), + v £E,- Avg (39) 


*Two parallel vectors do not necessarily have the same sense. 
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where z and (rv), describe conditions prescribed at the inlet and 
are constant on each streamsurface. Using Equation (36) this be- 
comes at once: 
[(V-i,wr)tQ=WxQ = РЕ, - у (ги), -Avp (40) 
If inlet conditions are uniform у 5. аай у (гу), are zero and all of 


1 


the vorticity is generated by the blades. In this case Equation (40) 





states that the vorticity vector, as well as the relative velocity уес- 
tor, lies in the blade surface 8 = constant; hence the term bound 
vorticity. This is the simplest example of bound vorticity. 

Forming the scalar product of the vorticity N and Equation 
(40), the left side is zero: 

Ave = MN: $ vLE- wtrv)] (41) 
Again this shows that the vorticity vector lies in the blade surface 
if the upstream conditions are uniform. 


Consider last the flow downstream of a system of blades. 


The equation of motion is simply Equation (40) with A = 0: 

WaN = РЕ,- v (гу), (42) 
If the conditions far upstream are uniform, v e and я (rv), are 
zero, and the vorticity vector is parallel to the relative velocity 
vector. Thus vorticity generated by a moving system of blades will 
be parallel to the velocity measured relative to the moving systern, 


regardiess of the shape of the blades, and regardless of the change 


in energy effected by the blades. This is merely the "upstream" case 





к 
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again if relative coordinates are used. 

The vorticity generated by a system of blades has certain 
geometrical prcperties. In the region of the blades in which the 
vorticity is generated the vorticity is 'bound' and lies in the binde 
suriaces. Downstream of the blades the vorticity generated by the 


blades is 'free' anc is parallel to the relative velocity. 
р у 





E) The Transport of Vorticity 

if a fluid particle has once acquired a rotation it tends to 
maintain this rotation as it moves along the streamline. In other 
words, fun vorticity is transported along the streamline with tke 
velocity of the fluid. The magnitude and direction of the vorticity 
will vary as the fluid expands or contracts, as the velocity vector 
changes direction, and of course under the action of non-conserva- 
tive forces. ‘ith uniform inlet conditions the equation of motion 
is given by Equation (39) with VÀ. and v (гу), =0. Taking the curi 


of both sides, 


242. V.víls Поги з рахе8в- 75 р: (43) 


This shows how the time derivative of the vorticity vector depends 
on the vorticity itself as well as on the velocity and the force field. 
If the fluid is incompressible the lust term of ¿quation (43) is zero. 
Because of axial syrnmetry, only the tangential vorticity is 
associated with the radial and axial velocities, while the radial and 


axial vorticity components are associated with oniy the tangential 








9 
velocity. The tangential vorticity may be regarded as an infinite 
number of vortex rings centered on the axis. The circulation /' 
around each ring is given by » multiplied by the cross-sectional 
area of the ring. As the radius of the ring increases the cross- 
sectional area decreases in such a way that the mass oí the fluid 
in the ring remains constant. It follows that for a given ring, /’ 
is proportional to D as the ring deforms. The tangential com- 
ponent of Equation (43) expresses the law governing the time rate 
change of ‘circulation connected wita a certain mass of fluid en- 


closed inanannular vortex ring: 





=/ /^ o v? Е ?//29 2 
"eV X = PIDE aL мог? 22:432] (44) 
From: the continuity equation, :zquation (21) 
/ 9 2 jel pp d ! 43 
— > rU + E&W = ۹م = کک‎ о Р 
r òr EIU P dr ie (957 
With this, iZquation (44) can be written: 
ou 2 v2 2 2742 44 
TE o2 [odiis м = 


The first term on the right is the rate of change in the axial direc- 
tion of the centrifugal force, which acts radially; the second term 
is the radial rate of change of the axial blade force; and the third 
term is the rate of change of the radial blade force in the axial di- 
rection. A little thought will show that these terms, in each case, 
represent moments tending to cause rotation of a fluid particle 


about a tangential axis. It is in this way that non-conservative forces 
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tend to effect a change in the circulation around a vortex ring. If 
the centrifugal force is conservative the first term is zero, and if 
the blade force is conservative A is constant, and the last two 
terms cancel. Then the circulation about the vortex ring is ccn- 
stant. 

Vorticity may therefore be considered as a property of the 
fluid which continually changes character as the fluid moves along 
the streamsurfaces. This concept of the transport of vorticity by 
the fluid particles gives physical significance to the non-linear 
interaction between the vorticity and the velocity and will be useful 


in understanding the various steps of the iteration process to be 





developed later. 
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V. Гн MATHEMATICAL FORMULATION OF 
THe IDEALIZED PROBLEM 

A) The Squations, the Variables, and the Boundary Condition 

Available for the idealized problem are six equations: the 
three equations of motion, Equations (20abc); the continuity equation, 
Equation (21); the isentropy condition, <quation (18); and the geo- 
metrical relation that the relative velocity is tangent to the blade 
surface, Equation (25). These six equations relate seven depen- 
dent variables, û, v, w, p, 0,À , £,--each of which is a function 
of the independent variables r and z. Clearly, since there are only 
six equations, one of these seven must be prescribed. The angular 
velocity of the rotor, w , is a constant parameter. 

in the direct problem the blade surface function f(r,z) is 
prescribed and there are six equations for the remaining six de- 
pendent variables. The necessary and sufficient boundary conditions 
are the sarne as for flow through the same channel with no blades act- 
ing except that a sort of Kutta condition must be applied at the trail- 
ing surface of the blade region. The boundary conditions will be 
discussed in more detail after the number of dependent variables has 
been reduced. 

In an indirect problem where the local rate of energy input, 
A (r,z), is prescribed there are again six equations and six depen- 


dent variables. The question of blade existence, so important in the 








a: 
indirect probiern with force components prescribed, does not arise 
here, for with the force field expressed in terms of the two functions, 
à and f, blade existence is assured. The boundary conditions, how- 
ever, are less understood than in the direct problem, particularly 
the conditions for f(r,2). lt is to be noted also that while the direct 
problem, with f prescribed, can be formulated as a problem in the 
meridional plane independent of A , the indirect problem, with A 
prescribed, is not independent of f. The two problems are therefore 
quite distinct. 

The direct problem is formulated beiow in terms of a stream- 
function for the meridional velocity. The discussion of characteris- 
tics and general properties in Part Vi also refers to the direct prob- 
lem. 

B) The Differential Squations for the Streamfunction 
The continuity equation is identically satisfied by the stream- 


function 4 i defined * by: 


EX cus o0 (47) 
ж P r az 
Е Е eoi ov (48) 
Es =F F Or 


, : XR XA. : : С 
For cases in which the fluid is incompressible ‘аке а, = 1,5 = |. 
The equation of motion is 


Vx N =vE-Avß (38) 


= In this definition a* and p* are considered as constants which de- 
scribe the initial total energy. H the initial total energy distribution 


is non-uniform, a* and p* are constants which describe the total 
energy at some representative point on the upstream boundary. 











у 
Upstrearn or downstream of the blades the force magnitude A is 
zero and since the total energy remains constant on each stream- 
surface 
2 = ور‎ (49) 


Written in intrinsic coordinates Equation (38) becomes 


a - | 
o v | = 225 (30) 
О» Л 425 


Using the definition of Op, anií2s , Squation (3), the normal 


component is 





T Ə rry О Е 2 
(Is =F on ~ dA (31) 
The other components show that ore is zero. 


“ince the total energy is constant on each streamsurface , 
Equation (34) can be replaced by 

ае. шу h-t AOE (52) 
V =- ln Jò (rv) ЖАУ бл» . 


and Equation (31) can be written for the upstream and downstream 


regions respectively as: 


а 7 9trr), DE, = 
PP EL ET 63 

M 2) 2(rv)e 2(rv, | 2E, | 
ИИС оз - TE (54) 


In the region of the blades A is not zero and the total energy is 
not constant on cach streamsurface. An expression for the tangential 


vorticity can be obtained from Equation (41): 


— 


N :v8 = De, 7 zz 5 - etr] (41) 


mxpressing the left side in cylindrical coordinates and the right side 
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n intrinsic coordinates there results: 














E Ee ET UT gabe etre] (55) 
But 14 2-47 5$ , Equation (3); and A7 g, $^ , Eq. (22b) 
so that 

а Ез ша ыш (56) 

A Ys ги ras 
and Equation (55) becomes 
IM - o rv E 
= r 2) ai zt- e rv] a 


Equation (57) expresses the tangential vorticity A in the region of 
the blades. It is significant that this relation is independent of the 
blade force magnitude A . 

All of the three basic equations, Equations (5 3), (54), and (57), 


are non-linear if the upstream conditions, E, and (гу), аге поп- 


1 
uniform, or if the fluid is compressible. If the fluid is incompres- 
sible and 5. апа (rv) are constant, the equation upstream is linear. 
It is simply n= 0, but the downstream equation, Equation (54), is 
still non-linear. However, the equation in the region of the blade, 
Equation (57), actually is linear if the fluid is incompressible. For, 
from Zquation (25) 


ку = ФУТ аг (би ғи) (58) 


so that in the region of the blades 


n = (kt - fef r0? r* (fou + tw) 
- 4 5 [&-wer)] (59) 





y= 


| the operator (Еа BE ды! represents differentiation 




















along the trace of a blade in the meridional plane. Equation (59) 

is clearly a linear equation for u and w when the last term is zero 
and f(r,z) is prescribed. 

in is written in terms of the streamfunction in &guations 
(5 3), (54), and (59), the three differential equations for the strearn- 
function are obtained. By direct substitution, using Equations (47) 








_ ¿ou aw 
ay SZ ЭР 
e 25 0*9 ER TEA (50) 
Pr or? F or 222 2 Рог год e oz 


The last two terms can be written rnore simply after a change of 


M independent variable as 
oy 1 2 Ф І 
ov P ò ðP _ af 1 OP (61) 


Or P Or 22 ед ^ on P дг? 


The derivatives taken normal to the streamsurfaces of 

` functions which are constant on the streamsurfaces can be written 
їп а very useful form by considering the streamfunction Y as an 
independent variable. For instance, since rv is constant on each 
streamsurface downstream of the blades, (rv), depends only on the 


streamfunction (p . Therefore 


O(re)u | e(rvk 99 ~- Alrvk/ з Р, (62) 
on dy on AY а, Ғ, 





and similarly for any function which is constant on the streamsuríaces. 


Using the relations developed in Equations (61) and (62) and 
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the definition of ф tne differential equations for the streamfunction 
can now be written for the three regions of the flow. Equations (53), 
(54), and (59) then become: 
Upstream of the blades: 


рн pes ЖИ а > kr у» (1 fln + (êr tg 25 u$s 8] 63) 


In the region of the blades: 
DAD Ba + (1+ RF) zz BLA a -(Z2R,-P Re) by 
= 227.27 + N2, - R2 (nF +10108). Ю2 фі, Un Pk 
22 2 о / г 
2 „г А er] (64) 


o v Ay? Cl y 


where R = гі, IE rí,, and the subscripts r and z dencte partial 





differentiation. 


Downstream of the blades: 
Prr B я Y, + wee 2 Фф» [In Pln -(2 h = зе" є, 
e 
The equations are rewritten below for the flow of an incompressible 
fluid with uniform inlet conditions in order to discuss the rotation- 
ality effects separately from the compressibility effects. For the 


incompressible fo Bid, a* and = are taken as unity. 


Upstrearn of the blades: 
Prr й = v. * 4; 70 (66) 
In the region of the blades: 


(/*Z*) - 27 фр: !! + RF] Vsz 
+ [ZE -R22 F] Yr- (ZR -RRA ж 2 rZ (67) 








Downstream of the blades: 

V, - Фа * zz = - [(rile- or?]Z (rele 168) 

With uniform inlet conditions the flow upstream of the blades 

is irrotational. ¿quation (66), which is clearly a linear partial dif- 
ferential equation, states simply that the tangential vorticity n is 
zero. Inthe region of the blades and downstream of the blades the 
flow is rotational. However, in the region of the blades the fluid 
is constrained to moveraent on the blade surfaces, and the vorticity 


15 'bound! and lies in the blade surfaces. Because of this constraint 





the equation for the streamfunction in the region of the blades, Equa- 
tion (67), is linear when the blade shape functions, R and Z, are pre- 
scribed. The equation for the downstream region, Equation (68), is 
a non-linear equation for | because y appears as an independent 
variable on the right side. However if the right side is a known 


function of position the equation is linear, although non-homogeneous, 





and can be handled in the same way as Equations (66) and (67). This 
suggests an iteration process in which the right side is estimated and 

` the equations are solved. A new right side is calculated from this 
solution and the process repeated until the desired accuracy is ob- 
tained. Actually (rv), is a function of ф only and the non-homogeneous 
part, the right side of Equation (68), can be evaluated as a function of 
Y and r for all the downstream region once the flow conditions at the 


iade trailing edge are determined. There are several reasons why 
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this iteration process is a practical and relatively simple method 
of solving this problem. ts part of Equation 
(68) "n a simple physical meaning. The right side is simply rn, 
the moment of the tangential vorticity about the axis, and is pro- 
portional to rer. It was shown that under the action of a conserva- 
tive force field the circulation /’ rernains constant along the stream- 
surfaces (Equation (46)). This suggests that a good first approxi- 
mation might be to assume that the circulation is constant on the 
streamsurfaces of an irrotational flow through the same channel, 
thus neglecting the interaction of n and the other vorticity compon- 
ents. actually a much better approximation can be made by evalu- 
ating the function (rv), as a function of y at the blade trailing edge 
using the blade geometry and the meridional velocity of an irrota- 
tional flow through the same channel. Onee (гу), is a known func- 
tion of y the right side is a known function of y and r. Here the 
interaction of the vorticity is accounted for, but the blade trailing 
edge conditions, which define rv(), are of course approximate. 
In the nurnerical method of finite differences used in the following 
examples it is only necessary to solve completely the final step of 
the iteration — . iach step only need be carried far en -gh 
to assure an improvement of the estimate of the right side. The 
application is relatively easy when the relaxation technique is used 


since the exact status of the solution is clear at all times. 
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When the fluid is compressible all of the equations are 
non-linear. The term on the right side, Yn (In eln , represents 
the entire eifect of compressibility if the flow is irrotational. For 
rotational flow another eftect of compressibility is to modify the 
effects of the rotation. ‘hus the density ratio (S yis a multiplier 
of the rotationality terms in Equations (63), (64), and (65). Again 
an iteration process is used, in which the term ty (ln Fl ry is esti- 
mated and a solution is obtained. This elution is then the basis 
for a new estimate of the density term, and the process is repested 
until the desired accuracy is obtained. , 

The conditions under wnich the above iteration processes 


are convergent have not been rigorously established. For flow in 





which a certain governing velocity is supersonic convergence is 
questionable. If stagnation points occur in the fluid the question of 
convergence is connected with the question of proper boundary con- г 
ditions and the fact that the streamsurfaces are characteristic sur- 
faces of the flow. Convergence will be discussed separately for 
each example solution. 

The equation for the flow in the region of the blades is 
somewhat simpler than shown in Equation (64) if radial blades are 
prescribed. Radial blades, blades which are generated by radial 


lines passing through the axis of rotation, are necessary in very 





-40- 
high speed rnachines because of the high centrifugal forces. The 
trace in the meridional plane of a radial blade is of course a radial 
line, hence f(r,z) is actually í(z) only. Then for radial blades: 
CEE Fee) GG 22), erc., (69) 
and Equation (64) becomes: 
ОТЫ аА та Yr + pzz = 22.5 ie 
* [1427] т,» talinela Hr [Ea 22) (то) 
The significance of the simplification is that the cross derivative, 
YP pz ‚ does not appear when the biades ere radial. 

The most general form of the differential equations for the 
streamfunction is that shown in Uquations (63), (64), and (65). These 
are the equations for the isentropic, rotational flow of a compres- 
sible fluid acted upon by an infinite number of arbitrary blades. 

The flow may be rotational and the energy distribution may be non- 
uniform at the upstream boundary. An energy relation which ex- 
presses density in terms of the gradient of the streamfunction and 

the total energy of the fluid must be considered with these if tne fiuid 

is compressible. 

C) The Isentropic Energy Öquation--Density as a Function of Mass Flow 

The energy equation, Zquation (32), can be written as 


+ и? с a = wry +£,- w(rv), (71) 





Let a, be the velocity of sound when the velocity V is zero. The 


constant a, will have a different value on each streamsurface if the 





"Ei 


energy on the upstream boundary is non-uniforrz;. .. non-uniform * à 


will be denoted by a prime, e.g. E ¡nen Lquation (71) becomes 





a? E / _ Си VES rV) 
a, 2 су. 2 


(72) 
or in terins of density 

/ 
Pl ш/ E шы ұ-/ wry - guía 2 
Р; = /1- 2 (араз) 2 SE Dé зема (73) 
Introducing the streamfunction defined by Equations (47) a 


(ae / ШҰҒЫЛ Гар a (О se) ee арен Р, 
е! Ра 2 p a= al 7 q* 2 


This equation gives the density in terms of the local mass flow, 
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pe tur git p? p? (75) 
and the effective tangential velocity 
2 
rer, (aa will) 
c = e (76) 


This last term includes the effect of changes of total energy as a re- 
sult of blade action as well as the actual effect of the tangential ve- 


% 


locity. 

“he solution of Equation (74) is oresented graphicaily as corn- 
putational curves, Figure éz. 

Two other useful forms of the energy eauätion can be derived 


{rom Equation (72): 
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D) The Matching of -olutions in Adjoining Regions 

[here are three distinct regions of the flow corresponding 
to the three differential equations, .2quations (63), (64), and (65), 
and solutions in adjoining regions must be matched on the connect- 
ing boundaries. [he streamfunction itself is always continuous 
across both the leading edge and the trailing edge boundaries. The 
Futta condition imposed at the blade treiling edpes to make the de 
lution unique requires that the velocity and pressure be continuous 
there. Consequently, the matching conditions at the trailing edges 
are that the streamiunction and its first derivatives, as wel! as the 
tangential velocity, are continuous across the trailing edges. This 
provides the means of evaluating the function (rv), for the down- 
Stream region. Since rv is continuous it may be evaluated just up- 
stream of the trailing edge by Equation (25), which can be written as 

ү <> ТОЕ еле ЖИ + kw] (79) 

the matching conditions at the leading edges are generally 
much more complicated. In the real problem with a finite number 
oi blades the fluid flows smoothly through the blade system. сл 
"stagnation point" occurs at some point near the leading edge and 
the fluid flows smoothly off the trailing edge. The flow is similar 
to the “low about an airfoil. However, in the idealized problem the 


number oí blades is infinite, and the velocity at the leading edge 
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vill generally be discontinuous. The flow is very similar to the flow 



















through an infinitely closely spaced lattice of airfoils, and the lead- 
ing edge discontinuity is the same as that occurring at the leading 
dge of the laitice. In problems of practical interest the pressure 
and tangential velocity will generally be discontinuous at the leading 
edge. Continuity of the first and higher derivatives of the stream- 
function depends on the shane of the blade near the leading edge. it 
can be shown that the first derivativesoí the streamfunction are cón- 
tinuous across the leading edges provided only that the leading edges 


lie in meridional pianes (Figure 4). 
К 






r n / y 
Leading 8 E 
Edge 
u 
AV 
Digure 4 
lhe Flow Across the Blade Leading 
Zdge 


It ia supposed that the tangential velocity v jumps discontinuously 
across the leading edge. This implies that a discontinuity in density, 


pressure, and velocity may also occur there. The first derivatives 





of the streamfunction are proportional to PJs > the local mass-velccity 


in the directions of the derivatives. па 18 of course continuous on 


both sides of the leading edge, and continuity of mass require that 
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Р) 19), = 12.) 19), in Figure 4. Cherefore pq,is continuous every- 
where, and the first derivatives of the streamfunction are continuous 
everywhere. sor blade surfaces which satisfy this condition, the 
matching conditions at the leading edge are that the streamfunction 
and its first derivatives are continuous there. The tangential, ra- 
dial, and axial velocity components and the pressure and density 
may be discontinuous. If the blades do not satisfy this condition 
the streamfunction derivatives will generally be discontinuous. 


2) The Boundary Conditions 


% typical region of flow is shown in Figure 5. The axially 


Downstream 


of 
Boun 7897 
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Edge 


Figure 5 
A Typical Region of Flow 


symmetrical surfaces AB and CD are boundaries of the channel and 
are streamlines of the flow. The upstream and downstream boun- 
daries are the surfaces ^C and BD respectively. The region of the 


blades or the body force field, over which the blade surface function 


f(r,z) is prescribed, is shaded. 
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* or each step of the iteration the equations for the stream- 
function, Equations (63), (64), (65), are non-homogenzous linear 
second order partia! differential equations. The necessary and suf- 
ficient boundary conditions for this problem are well known. ¡he 
streamíunction $ , ora linear cornbination of its derivatives, Y, 
and Y, , must be prescribed at every point on the boundary. The 
boundary conditions on tne streamfunction for a typical region of 
flow (Figure 5) are: 

Ш = constant on surface AB 

ф = constan: on surface CD 

DH = Y (r) on surface áC 

E = 0 on surface BD 
The functions У, р, апа р, which can be preserbed only once on each 
streamsurface, are prescribed at the inlet, on surface 3U. 

st the leading edge the matching concitions are that the 
strearnfunction and its first derivatives are continuous. At the 
trailing edge the streamfunction and its first derivatives, as well 
as the pressure, density and tanpential velocity, are continuous. 
¡his is the Kutta commisi. 

The iteration process is based on an elliptic partial differ- 
ential equation. in Part VI it is aiso shown that tke complete non- 
linear differential equation, which inciudes tus compressicility 


term, becomes "hyperbolic" when the meridionil velocity or, if 
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blades are present, the relative velocity, becomes supersonic. 


LA 


The boundary values for this case are quite complicated and solu- Y 


tions may not be unique. In all the examples solved here the gov- 


Р” d 


ж 
erning velocity is subsonic, although the total velocity may be 4 [ У | 
| р 

transonic, as in Part X. 
Е) The Limiting Flow Far Downstream of the Blades 
If at a finite distance downstream of the blades the channel 
boundaries become concentric cylindrical surfaces and if these 
surfaces extend unchanged downstream to infinity, general quali- 
tative statements can be made concerning the limiting flow at in- Ж 
- | 
finity. It is supposed that the total energy of the fluid is uniform ж 


and that the vorticity and tangential velocity are zero on the up- 
| - 





stream (inlet) boundary. The energy of the fluid is changed and 
vorticity is created by the action of an infinite system of blades 
rotating at a constant angular velocity W . Nothing is said about 
the channel boundaries upstream of the cylindrical region or about 
the shape of the blades. 

It is clear that in such a region the flow is independent of 
the axial coordinate and the radial velocity is zero. Then radial 
equilibrium of the — force and the pressure gradient 
force requires that (frorm Cquation (20a)} : 
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The vorticity, being generated by rotating blades, is parallel to 


the velocity relative to the blades. Therefore from Uquation (42): 





E IV e ^v s (81) 


From quations (80) and (31) and the isentropy condition, Equation 


(19), the following relations are derived: 











Ze. = 
Eee T (32а) 
da? Е 
А т (Y-1) -Е- (82b) 
2 
ect. _ bs ағ er sex) 
Fal a = 
dM? _ 2 = 1-1 142 Dr dry 
dr °° a? eg ie) (82d) 
The energy equation, ¿Equation (72), can be written as 
2 
0-1 e “о -/ wr 
Using this, 
am? E E 3 Ll zur 
adr Ба Ко су 2 (82е) 





If rv is a known function of r or of y in this region the flow 
can be completely determined by the above equations. It can be 
seen directly from <quations (80), (82a), and (82b) that the pressure, 
tne density, and the local speed of sound are minimum on the inner 


boundary and increase monotonically to a maximum on the outer 


Ше 




















boundary. rhis general result is completely inde pendent oí the blade 


LE 
bed bu 


shape ¿nd the channel shape. lí the blades are stationar y it follows 
from Zquations (82c) and (82e) that the total velocity and the corres- 
ponding Waci number are maximum on the inner boundary and de- 


crease monotonically to a minirnum on the outer boundary. The axial к. 


Ц 

velocity may have a maximum or minimum anywhere, depending on + We 
y 

, 





rv, but if rv is monotonic then the axia! velocity is also monotonic. ) і 


Іп order to examine the total mass flow the channel is given | 
the sarne cross-sectional area far upstream and far downstrearn. 
The total energy of the fluid is constant. If the tangential velocity 
is zero upstream and constant downstream the total mass flow egua- 


tion is 


ср ае 





2) = Га Wa rdr (83) 


where ( ) and ( ); denote values far upstream and far downstream 
respectively. 
Using the isentropic ge ton, 

2 = (1 - Е: ER Ж” 7 


and integrals of Fquation (82), Zquation (83) becomes 


f / 
| / = 2 yey. 2 2 2 2 
E үл М" J” М «ef "y? см, - Enri) Jut Seine rar (84) 


where Mf s E] + Ж is the maxirnum which occurs at the inner 
2i 





radius, taken here as unity. 


«he solution of '¿quatiocn (84) is shown in «igure 43 where 
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P is plotted as a function of /V for several values of A : 
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These results are expressed in terms of local Mach numbers in сірге 


2 
44, where the axial Mach number (Y) , is plotted as a function of 





2; 
2 
the inlet hiach number м,? for several values of == . Itis seen 
2; 


that for a given inlet Mach number, M, 2 , and a given deflection, c і 


a; 
there may be two solutions for the maximum axial Mach number 
- . v 2 е 8 в. е Й 
downstream, but if 2) is above à certain critical value there will 
2; 

be no solution. When two solutions exist one value of the axial ve- 
‘locity may be subsonic and the other supersonic, or both values 
5 Y ғ ` 1 T 4 , Ww. vv 
may be supersonic. It must be understood that (2) ала 

i e 2-/ As m y 
are maximum values on the downstream boundary and that when 
these are sunersonic the minimum velocities, and, even the mean 
velocities, may still be subsonic. 


These results depend only on the limiting flows far upstrearn 


and far downstream and are independent of the intermediate flow. 





T 8. 


It is assumed, however, that the intermediate channel and the blades 
are such that an isentropic transition from the upstream flow to th 
downstream flow is possible. Nevertheless, this analysis shows 
that for a given inlet Mach number there is à maximum tangential 
velocity that may be generateu by stationary blades. No isentropic 
solution exists for tangential velocities exceeding this maximum.. 


Similar results would be obtained for distributions of tangential 


velocity other than constant, and for rotating blades. 
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VI. PROPI2TIES OF THZ PRANSONIC FLOW 
Or аА PERFECS FLUID UNDER THE ACTION 
OF 42 POCUDO-CONSERVATIVE FORCE FIELD 
А) “he Characteristics of the Problem--The Governing Velocity 
ihe applicable equations upstream and downstream of the 
force field are shown in Equations (19), (20), and (21) with A =0. 
Adding to these the four differentials for u, rv, w,ln p, and using 
Equation (19) to eliminate the pressure terms, eight algebraic equa- 


tions for the eight partial derivatives are obtained: 


ги да еле oy Ove 
“сз; ^ M Ty а = = 
дг г, = 
or d o2 a 

ду дум ты оГ М c 
A E CATRE. © 
Ju ow / gf ‚ дё _ A 
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dr OU, ez OU = AU 
Zr деу, ez orv 22 ү гү 

д д 
dr èw + حلب‎ 20 = dw 

О 2. 
f l / ef = x 
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Ihe differential equations for the characteristics are obtained when 


the determinant of the cocfíicienis is set equal to zero: 











“¿xpanded, this equation 15: 


м?/м? ага“, 22 и -а ат ат +» [6utw*/uw+w*Ya?]artadz? 





- 2uw(2u* a" ]araz? + u?lu-a)dr* = 0 


which on factoring yields: 





d и |2 | 
о зо (87) 
dr | uwzayut+w?-a? 


 Equations (87) and (68) are the differential equations for the physical 
characteristics of the nroblem. These equations represent, respec- 
tively, the streamsurfaces and the 'meridional Mach surfaces". 

Thus the character of the flow upstream and downstream of the blades 
depends on the meridional velocity, not the total velocity. That 15, 
the problem is "eliiptic' with respect to the variables u and w if the 


meridional velocity is subsonic, and is "hyperbolic" if it is supersonic. 














The veiocity which determines the character of the problem is 
called the "governing velocity", thus the governing velocity up- 
stream and downstream of the blades is the meridional veiocity. 
in the region of the blades the applicable equations are 
those given above but with A# 0. In addition ¿quation (25), which 
states that the relative velocity is tangent to the blade surfaces, 
must be included. Zquation (19) can again be used to eliminate 
the pressure and Equation (20b) can be used to eliminate A 
Because of the constraint imposed by the blade surfaces (Zquation 
(25)) it is also possible to eliminate the angular momentum rv. 
The blade surface function, f(r,z) is prescribed. Then, including 
the differentials for u, w, and ln p, there are six linear algebraic 


equations for the six partial derivatives: 


ES two - ара + wz | +a? 
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du 2м ‚ ә? | oP И 
CAS E И ааг >= م‎ 
и (89) 
dr az 2% =н 
r 
с/г ож + Az = dw 
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setting the aeteri::ingant of the coefllicienfsequal to sera tne 


following equation 15 obléined: 


wj (I IIA ng) +) ر‎ Із( ӨЗ uw? (вФіа-2РФМДа аға: 


2 [o 9^ 27) aw -[ IH Z)w-2RZ uJa?}ar dz * 


= u {+ R*+2 u [(ı+22)a*] dz = O (90) 
On factoring this yields: 
Wdadr-udz =o , ame 
/##'42* py з HR?) dr? - B (8 2^ uw + LES | RE Zid raz 
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she differential equations ter the physical coaractoriztics are 


A (91) 
Cu = (ек uw +RzZatzra/(mr2 ЫЗ (Rarzw)? a] (42) 
dz (i+ R24.27) #- (кю?) а? | 


where Mu+Zw is the relive tangential velocity, (v-w kr). 

ihe characterigtics o1 the flow in the region of the bleues 
are the streamsurfaces (lquation 91) опа the surfaces given vy “qua- 
tion (92). ‘/heth@r or not the charicteristics Given by Equation (92) 


are real depends on the total relvtive velocity, not on the me: *idional 
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component cf the velocity as in the upstreatn aod dcwnsireuirm regions. 


Che governing velocity in the region of the blaúes is therefore the 
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velocity relative tu the blades. 
Phe form of cquetion (92) suggests the following transform- 


ation: 


//»02,22 //+# @®*=2* E, 


а = u ; w= 
"TE / ur 
(93) 
Ar = ъв? ат ; AE sy ES oz 
After this transformation Equations (91) and (92) become: 
ez. u 
ae W 
a= 222 с? 2.2 DRZIMR%Z _ _ IRZ _2 {са 
зада + UL A a? (91) 
E Р Пн ИН etz) (ена) 
ez. E 


This is no great simplification except when R or Z is zero everywhere. 


Por radial blade R is zero, the transformation is: 


4 = а : Ж т/у/ғ22 уу 
(95) 
іг = іг; dz=VJ//+z2 az 


and the equations for thu carracteristics, in terms of the transformed 


variables are: 


ero uw 

сар (96) 
ФЕ üUw-rzta/a?,w?-a* 

dz М2-су2 


Tkus the characteristics for the transformed problem are the same 
as for two-dimensional flow with velocities U and W in the F, Z space. 
5) The Complete Non-linear Forms of the Equation 

With unifcrm inlet conditions the equation for the strean.- 


function for cornpressible flow downstream of a system of blades 
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rotating at angular velocity W is (from Tquation (65) ): 


Prr- F Pr н ае = 4 (10 0), [DO [Lt - 2] L, вм (97) 
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The density derivatives can be written using Squation (78) as 
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Then, eliminating the density derivatives from Equation (97) and 
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introducing the velocity components, the complete non-linear form 


of the equation for the downstream region is obtained. 
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The non-homogeneous part, "the right side', represents the effect 
of rotationality. The terms containing the derivatives of > are 
connected with the fact that part of the total energy occurs as kinetic 
energy due te the tangential velocity. it is clear that the effect of 
rotationality becomes zero if the meridional velocity is sonic, and 
that the rotationality effect on one side of the sonic line is exactly 
opposite the effect on the other side. 

The equation for the streamfunction for compressible flow 


in a region of radial blades is (fren: Equation (70) ): 
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The density derivatives, írom Equation (78), can be written as: 
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The complete non-linear equation for the region of the blades is ob- 
tained by eliminating the density derivatives from Equation (100) and 


introducing the velocity components in the coefficients: 


(1-5) (1427) pp + 242 11427) pa + [1-14 27128 J Yar 


dr z2)e (12) Les [24 (142) w lA 
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С) The Cushioning “ction between Subsonic and Supersonic Regions 
The compiete non-linear equations (2quations ($9) and (102) ) 
are obviously quite complicated. It is believed, however, that the 
transformation suggested (Equation (95)) would lead to some simpli- 
fication of quation (102). Unfortunately, time has not permitted 
a more thorough investigation of the possibilities of this transform- 
ation. 
Examination oí Squations (99) and (102) and the characteris- 
tic equations, Equations (88) and (92), reveals that the non-homo- 


geneous part, the right side of the equation, changes sign whenever 
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the character of the flow changes from "eiliptic" to "hyperbolic". 


Furthermore, for flows of practical interest the vorticity term 


[2 | у NONEM дк mme oue E 
Fe а. Ay AY I 2 о, Px 


+ 





will have the same sign throughout the flow region. The right side 
of the equation can be regarded as a forcing function acting to dis- 
place the streamlines. In regions where the flow is entirely sub- 

sonic or entirely supersonic the forcing function will have one sign 
throughout the region. Kut if the flow ts transonic, i.e. part sub- 
sonic and part supersonic, the forcing function has one sign in the 

Subsonic part and the opposite sign in the supersonic part. There- 


fore the deflection of the streamlines, brought about by vorticity, 







will be less if the governing velocity is transonic than if it is en- 
tirely subsonic cr entirely supersonic. Consequently the meridi- 
onal flow will be "smoother", for the vorticity in one region tends 
to counteract that of the other region. It is believed that this mu- 
tual cushioning efiect is the explanation of the phenomenal effic- 
iences observed in compressors in which the relative velocities 
in 2 region near the tip are supersonic. 

The fact that the forcing function, the non-homogeneous 
part of the equation, contains a factor of the form (-м,2 , Where 


ім 


g is the "governing' Mach number, leads to the following general 


conclusions: 
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igth of vorticity at a certain point in the flow is sero 





when the governing velocity is sonic at this point, has one 
sense when the velocity is subsonic, and the opposite sense 
when it is supersonic. 

2. the deflection of the streamlines brought about by vor- 
ticity in a region is less when the governing velocity is 
transonic in the region than if the velocity is entirely sub- 
sonic or entirely supersonic. 


These conclusions are confirmed in all of the examples of 













compressible flow which follow. In fact, it seems that even stronger 
statements should be made. A comparison of the deflection of the 
streamsurfaces brought about by vorticity indicates that the deflec- 
tions increase from zero to a maximum and then decrease as the 
vorticity (the strength of the actuator) is increased from zero. 


The governing velocity is subsonic in both examples. 
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vl. PHE FORMULA TIGN OF TH ги 
DUP ERENCE PROBLEM 


ТТЕ 


fhe finite difference problem differs from the differential 
equation probiem inu two ways: in the finite difference problem, 
the boundary conditions are prescribed at a finite number of points 
on the boundary and the desired function is to be determined at a 
finite number of points within the boundary; the differential egua- 
tion is essentially replaced by a number of simultaneous approxi- 
mate difference equations which are themselves only approximately 
solved. However, the set of points at which the desired function is 
sought can be rnade very dense, and the system of difference equa- 
tions can be solved as accurate!y as desired. Cherefore, barring 
unaccounted for singularities, any degree cf agreement betwee 


the two solutions may be obtained. The two problems are neverthe- 





less quite distinct. fhe finite difference problem is solved by "outh- 
well's relaxation method (Reference 8), which is defined as "a sys- 
 temätic sequence of localized changes of the wanted function that 
steadily brings the 'residuals' toward their desired value." 

ihe fact that the physical asnects of the problem are always 
evident and that the status of the solution is always apparent is a 
great advantage of the relaxation method. This is a great aid in 
solving the non-linear turbomachine problem, for the action of the 
vorticity and the compressibility is apparent as a force" which 


causes the streamlines to deforyn fom their incompressible, irro- 
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tational positions. 
Yhe power oi the relaxation method is indicated in the fol- 
lowing quotation from Southwell (Reference 5, р. 3): 
"(The) use of finite differences is no. new, nor (is the) 
evaluation of the wanted function at nodai points of a regu- 
ler net. bu: in concenirating attention on the data, amu in 
recognizing that these are never exact, they subordinate 
inathe matical to physical aspects in a way that can alter 
drastically the course ої а thecretical research. Dis- 
carding orthodox for relaxation methods, an investigator 
finds his outlook quite transformed: full scope remains 
ior ingenuity and special artifice, but any problem that 
can be formulated can be solved.’ 
It is not clear how solutions which are unstable can always be found 
by the method of finite differences. Intuitively it seers that if the 
configuration is unstable the disturbances inherent in the reiaxation 
process may be amplified by the relaxation, thus making converg- 
ence апеєіісласіе . 
A) The Difference ._quations 
ihe most general form of the differential equation for the 


streamfunction is that shown by Equation (64): 


Афр + 28 рг + Сф + 2 = О 
(193) 
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‚here 2,B,C are function 


Шр Шо т, апі 2. As written for each step of the iteration with 





; of r and z, and D is a function of 


the non-homogeneous part known, Equation (64) is an elliptic eaua- 


tion, for 


NB a eee 50 
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and can therefore always be transformed into the normal forrn 


фух ож Pyy «20 (х,у, Ф, Жі 5 О (104) 


ya 


When the cross derivative term /,; occurs, as in Zquation (64), the 
difference equation may be written for either the original equation 
or the transformed normal equation. Because of the added difficulty 
of matching solutions when one region has been transformed, the 
difference equations are written here for the original form. 
If the streamfunction Y , which is a solution to a partial 
differential equation such as Squation (64), is regular in some neigh- 
borhood of the point (rz, 2.) , it may x represented in this neigh- 
borhood by the Taylor's expansion: 
pinzi = P + 1-6] + dz (2-20! 

+37 [brririlie 2 (е) о) + pza (2-2o) 


mM. (105) 


where Y, is the value cf Y at(r,¿,2,) and ali the derivatives are 


-ъ - 


evaluated at (rı T . There- 
ry s + E 
fore if (ro. 2.) is one point cf 


& square lattice with an in- 





terval between points of ô, | 
Figure 6 


The Square _attice cr Net 
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the value of y at the surrounding points can be expressed by mtans 
of =quations (103). If fourth and higher powers of & are neglected 
the result is a system of simultaneous linear algebraic equations 


which can be solved for the derivatives to give: 


Y, -Va- P4 Pa е Шъ ج‎ ф, - фіз 
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On eliminating the first order derivatives it follows that: 
V, + bo + Ya +Ya-Abo= [Yrr+ zz) 5 
Ш, + db, + Y, +De-Ado = (Yrr+ Y22) 5 
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(107) 
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[he difference equation is obtained when Zquation (106) is 
substituted into the differential equation. Thus for Dquation (64), 
which is 

/ 
(122) 0,-2Р2 г (100) 0, + {22,-#®2, у) ф„-(2В-®ЮЮ,) „= @ (108) 
the difíerence equation is: 
5 
[irl PRIME] by + өг2йгр өлгі 610, 
¿(RAZR RRE] Ys + [127]- [22-22-21 ] dy 


"ERI. ode) – 41 518247 0, = QS? ^ — (05) 
where 


R= £ z+ [ltl $,-R2 Ye] (In A, + [RI by - ez, I(In d, 





£, / 
НЕЇ є 257. E 





-53. 
When the blade shape function f(r,z) is prescribed all the coeffic- 
ients can be determined. Equation (109) is then the difference equa- 
tion applicable in the region of the blades. The residual, A 5“, 15 
a complicated function of position which is estimated for each step 
of the iteration by means of the solution of the previous step. 
In the same way the difference equation upstream of the 
blades is, from Equation (63): 
pW + (lÊ) Va + Фа е а - 44о - @5° (119) 
where 
R = 4, (1те), 4E [77 z. E - ir) 
The difference equation for the downstream region is the 


same as Equation (110), but the residual, from Equation (65), is 


given by 
@ = e ДБ). AE E, E zz 
2/27 2d £ 2 а ан) 
(EJ [772 RE - mE ар a 


The difference equation for each of the three regions is 
shown in Equations (109), (110) , and (111). A general form might be 
written as 

Dily] = SR (112) 
Here D, (4) represents the leít side of one of the difference equa- 

th 


2 
tions written for the i` point of the net and R is the residual 


there. D,( ) is the "difference operator’, and is itself a function 











of position. fhet is, the coefficients of the left side of the differ- 
ence equations have different values at each net point and these 
vaiues depend on the prescribed Shape of the blades. 
if the blades are prescribed ag radial, the expression for 
D; is more simple than in the general form (Equation (109) ). For 
radial blades 
Did] - Y, "L9 22)- (1-22) E- Ta + p3 
+ [liz 01-2275] Ye - 4[1+327 b, 


in the region of the blades. The (A; are correspondingly more 


(113) 


simple for radial blades. However, the compiexity of @ is pri- 
marily due to rotationality and compressibility. For example, if 
the fiuid is incompressible and the inlet conditions are uniform the 


right side of Zquation (64) becomes simply 


A, = 2r 2Z (114) 


J 


regardless of the blade shape prescribed. 

-pecialized forms of the difference equations will appear in 
the various examples, Parts VIII, ix, and X. The general form of 
the difference operator is written as the sum: 

D, (¢) = 2 а, Vj (115) 
The subscript i denotes the point of the net for which the operation 
is performed md the subscript j denotes the values of the stream- 
function which enter into the difference equation. Thus for the 


simplest difference operator, from Equation (110): 
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D; (4) = Wy, + (1-8, ) by2 я sr (il ti з 4% 


о, 5 ар = / D Ao = Т й (1 16) 
ie E г == 8 
C2 = 2; 7 C, 7 /* 2r. 


The a, are called "influence coefficients" in that they indicate the 
/ 


2 1 ° « 
change in 8°; , corresponding to a unit change in фуу ; ог =; 


represents the change in фії due to a unit change in the "forcing 


function" S R, when all other Y are held constant. The A; from 
J 


Zquation (116) are shown for sorne point i: 
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The Nomenclature cf the Difference Cperator 


"igure 7 
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rhe difference equations must be satisfied at every point 
This is accomplished by southweli's relaxation tech- 


in ihe net. 


nique, whereby the streamfunction is modified until the desired 
Ф 


residuals are obtained. 


The soluticn of a finite difference problem is generally 
not the same as the solution of the corresponding differential equa- 


tion problem. Ihe amiouni of the discrepancy depends on the nature 
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of the function involved and оп the size of the net used. There- 
fore it ig only necessary to obtain an accuracy in the solution of 
the difference problem commensurate with the discrepancy already 
involved. 
E) The Iteration Process 

it is not practical to attempt to formulate a generai set of 
rules which will govern the iteration process in all conecivable 
problems. It is, however, appropriate to discuss generally the 
two separate effects--compressibility and vorticity. For the initial 
step of the iteration it is necessary to estimate the right sides of 
tie equations. Consider, for example, the right side of Equation 


(65): 








2 
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‘he manner in which rv(y) may be estimated has already 
been discussed. in the examples solved the first estimates of 
ку(ф), со obtained, were accurate to within three percent of the 
final values (Figures 2l and 23). Ihe complete effect of rotation- 
ality, represented by the term in the bracket above, can likewise 


be expressed as a function of Y aná г with comparable accuracy. 


In the relaxation precess it is possibie to improve the desired 
residual each time tae strear..function is improved, so that the 


relaxation and the iteration are actually perforrned simuitaneously. 


EM. 
Periodically rv(Y ) must be recalculated and a new, riore accurate 
function of Y aná r obtained. 
rhe compressibility effect cannot be handled so easily. If 
the Mach number is low and the rotationality effect is strong, аз іл 
the examples of Fart IX, the first term, di, (In e), is much smaller 


than the last term and can be disregarded for the first estimation. 


In this case the effect of compressibility is not completely neglected, 





but appears in the (| factor multiplying the rotationality term. 
it is true that under these circumstances the main effect of Сотпргез- 
sibility is its "influence" on the гоба спа гу. Гог large subsonic 
Mack numbers this approach is not possible. The first term may 
De of tie same order as the second. In fact, in the examples of Fart 
%, where a constant tangential velocity was prescribed at the down- 
streara boundary, the two terms were exactly equal at the point on 
the downstream boundary where the axial velocity became sonic. 
“his can be seen by writing ¿quation (65), with ш 2 0 and Jz4 - 0, 
as 

(Е) е (Рм) < 0 
it is well known that if the energy is uniform the local mass flow 
із а maximum when w is sonic, һепсе is zero when w is sonic. 


inother viewpoint is to regard the residuals as forcing functions 


which act to displace the streamlines. Since the axial velocity can 








Ое 
only be sonic at à 'inroat", the forces must act toward the sonic 
line so as to contract the streamlines there. Ко residuals, having 
a different sense cn each side of the sonic line and being 'continu-- 
o0us'', must therefore be zero on the sonic line. ‘shen vis act con- 
stant or when the energy is non-uniform, Pw will not necessarily 
be maximum when wis sonic. The exact circumstances under 
which the residuals change sign were derived in Pert Vi, where it 
was shown trat the residuals are zero on the line separating the 
"hyperbolic" and the "elliptic' regions, or, in other words, where 
the ‘governing velocity’ is sonic. 

It is extremely irnportant to use every means available to 
make the first estimate of the compressibility effect as good as 
possible. The rapidity of convergence depends on the judgement 
used in the first esticnation and each successive approximation. No 
general rules can be stated. In regions where the meridional velocity 
is nearly sonic, the streamsurfaces obtained during the iteration may 
be such that the channel between adjacent surfaces is ''choked'’. 
When tnis happens the rnass flow Рад, falls to the rignt of the density 
curve in Uigure 42 and the density is imaginary. If it is not possible 
to Edjust the streamlines so as to have real densities, then it may be 
that the governing velocity is supersonic and the ' elliptic”" forrn of 
the equation is not applicable, or even that no isentropic solution 


exists. if there is reason to believe isentropic solutions exist, 
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the "complete non-linear" equations derived in the appendix may 
oe used for the transonic problem. Convergence is always slow 
when the governing velocity is nearly sonic, for the p vs pegcurves 


have nearly vertical tangents then. 











Vili. PINEXANPILE OF INCOMERE£TIML E FLOW IN A MIXED 
Im СОАК ТАОМ WITE BONS SNAPE PRESGRIBEL. 
This analysis of the flow in à mixed-flow comnressor with 
blade shape prescribed is presented as an exarnpie of the three- 
dimensional moticn of an incompressible fluid under the action of 
a systern of rotating blades. first a soluticn is cbtained for the 
irrotational flow through the channel with no blades present. This 
serves as @ basis of comparison for determining the additional! 
velocities induced by the blades, and is alsc useful in making the 
firit approximation to the rotational flow when blades are present. 
the particular channel chosen is shown in sigure 8. The 
curved part of the channel boundary can be expressed analytically 


as 


ғ 
22. (117) 


7- 
Pa“, 


= 
po 
where ri is the initial radius ‚г is the final radius, end L is the 
length of the curved portion. The coordinates of the channel boun- 
Garies are given in Tabie I. 

In this example the total energy is uniform and the vorticity 
is zero at the inlet, station 0. More precisely, the tangential ve- 
locity is zero, the meridional velocity is axial, and the pressure 
is constant on the upstream boundary. In the blades and downstream 


of the biades the total energy is non-uniform and the flow is rota- 


tion:!. The downstrearn boundary conditions are essentially those 
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that might be expected if the chennei extended downstreara to int 
iinity: the flov is independent of the axial coordinate. Thus the 
radial velocity is zero and the pressure gradient forces and the 
centriiugal forces are in balance on the downstream boundary. 


A) Irrotational Flow with no Blades Present 





For irrotational flow with no blades present the governing 
differential equation for the streamfunction (Equation (63)) is simply: 

Prr -+ Pr + رل‎ = О (nn 
¿nd from Equation (110) the corresponding difference equation is: 

b+ 11-6) da + (14) - 44, = O (119) 
fhe boundary conditions are: (Y = constant on the hub and shroud 
contours ,-24 =w=/ cn the upstrearn boundary, and 0, = и = О 
on the downstream boundary. Upstream y, is taken along the boun- 
ary, so that Y czan be determined by simple integration and pre- 
scribed directiy. Downstream the normal derivative is pas wer bed 
and Û cannot be determined by integration along the boundary. The 
boundary value problem is therefore of the mixed type. However, 
Since the flow far downstream is entirely independent of the axial 
cocrdinate all axial derivatives vanish and the differential equation 
(Equation 118) can be integrated easily. By this means y itself can 
be prescribed over 3il of the boundary. The mathematical problem 
for the irrotational ilow is therefore relatively simple and consists 


сі satisfying simultaneously th» difference equations for all points 
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Le net. 

The flow region is divided into a net or lattice and influence 
coefficients are calculated for each point. The net, the influence 
coefficients, and the boundary conditions are shown in Figure 9. 
Che influence coefficients near the boundaries were calculated in 
the same way as Zquations (106) and (108), but with & different 9 
for the "short legs" of the net. 

For a first approximation to Y the streamsurfaces were 
assumed to have the same shape as the boundaries. This is equiva- 
lent to assuming that the axial velocity is constant on each radial 
line. The actual residuals g for the smooth values of 2 меге 
calculated by means of the difference operator with z high degree 
of accuracy. These served as basic values of W and Rq . The 
change in Y to make the residuals as small as possible was then 
Getermined by the relaxation process. As a final check residuals 
"were again calculated by the difference operator. 

The final values cf the streamfunction and the corresponding 
actual residuals are shownin Figure 10. Upstream of station 12 
the smallest residual may be as large as 2 and downstream as large 
as 5. This error corresponds to an error in the solution of the dif- 
ference problem of about cne part in 2000. 

The axial and radial velocities for several representative 


stations are shown in Figures ll and 12, respectively. fhe abscissa 








E. 


of Figure llis the deviation from the mean axial velocity: 


W- Wm 


aw = 
Wo 


where w, is the inlet velocity and w_ is the mean axial velocity. 


The ordinate is the non-dimensicnal radius: 


т. E f- Гу 
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where r, is the radius of the inner boundary end rg is the radius 
of the outer boundary. 
The abscissa of Figure l2 is simply: 


ЕЗ 
u~ We 


Lines of equal pressure are shown in Figure 13 for flow 
through the channel with no blades present. The decrease in pres- 
sure occurring far downstream is due to the decrease of the cross- 
sectional area of the channel and the corresponding increase of the 
mean velocity. This should be taken into account when evaluating 


the pressure increase when the rotating blades are present. 


B) Rotational Flow with Blades Present 





notational flow through the channel of the previous example 
is considered here. The vorticity is generated by a rotating sys- 
tem of blades which acts on the fluid between stations 7 and 14, Ra- 
dial blades are prescribed, hence the blade sbaps function f(r ,z) 
depends only on the axial coordinate z. Разед оп Ше irrotational 


meridional velocities of the previous example the function f(z) was 













the rate of change of 
the moment of angular rnomentunmi rv, was very sinall near the 
leading and trailing edges. The discontinuity in v across the lead- 
ing edges cannot generally be avoided, however the blades were 
selected such that the jump in v is positive. The rate of energy 
input is indicated by the slopes of the curves in Figure 20, where 
the desired and the actual > on the radius r = 2.5 are cornpared. 
The blade shape function is chosen “o give the desired um on the 
basis of the irrotational meridional velocities. Ihe actual A is 
determined from the final meridional velocities with the prescribed 
blades present. The shape oí the blades is shown in Figure 8 and 
the blade surface function is given in fable i. 

The applicable differential equations for the three regions 
cf flow were derived in Part V-E. 


фур E p O upstream of blades (66) 


E. - (1-22) pr + Pez 


а Ра Є in region of blades (27) 
Ur E №, ЕЕ m - pra" downstream of blades (68) 


where rv = rv( Y) only and is evaluated at the trailing edge by Equa- 


tion (56), which nere becomes: 


[rvh = wrt [7 +( Z), | (120) 














The corresponding difference equations in regions where the net 
interval is 0 were given by Equations (110), (11i), (113) апа (114), 


which for the three regions become, respectively: 





V, E dis + (1+) by - 40,0 (121) 


b+ [-E)+ (+2127 by + by + [VB )+(1-£)27 be 
-2[1+222] 4, = 2.02 52 (122) 


Фе В, + Ys + (+ Fl da 44, (rv- or] ST" s* (123) 
where rv(%) is evaluated by Equation (120). 

Lhe influence coefficients for chese equations and for the 
equations which hold near the boundary where $ is not constant are 
shown in Figure l4. 

The boundary value problem is the same às in Part VIII-.., 
except that here the streamfunction can nct be prescribed exactly 
on the downstream boundary. Instead new downstream boundary 
conditions must be found for each approximation to rv(y). This is 
easily ure by one-dimensional relaxation on the downstream boun- 
dary. The first approximation to the downstream boundary values 
was based on rv determined from the irrotational flow. For com- 
parison, the initial and final downstream boundary values are: 

r 3 29 3.0 ст» 279 2599 2.0 
initial Y 0 973 1557 2654 3367 4000 


Pinal Ф 0 977 1863 2660 3371 4000 








ва 
“he parameters ct the problem, based on tbe сплохітога 


radius oí 3.425, are: 


Pressure coefficient: т 2 3655 
Flow coefficient: = = 0.769 
Root=tip ratio: p^ = 3 upstream 
= 1.625 downstream 
f = 2 upstream 


п 


‚625 downstream 


Blade aspect ratio: 0.232 


static pressure increase: (22) - (ЕЕ z 1.570 at tip, 
blades 


Фо Fo} no blades 





= .756 at root 


three-dimensionai properties of the flow, and blades of very low 
aspect ratio were nrescribed in order to stress the importance of 
blade geometry. 

there was nothing unusual in the analysis. Relaxation and 
iteration were performed simultaneously, and tho function rv(4 ) 
oniy had to be determined once after the first estimation. the 
initial and final values of rv(¢) are plotted in Figure 2l. The tinal 
values of the streamfunction and the corresponding desired and ace 
tual residuals are shown in Figure 15. The axial and radial veloci- 
ties are plotted fer several stations in Figures lo and 17, using the 


same variabies as already described. iines of equal pressure are 


= 765 
sown itu igure 19 апа lines of equal total energy in Figure 19. 
The pressure was calculated by a relation derived from Zquation 


(31): 
DES. гу | Po u* sv» uw? 
z Pw? EE Мо 2 (124) 


A comparison oí the two solutions, with and without blades, 
does not lead to any extraordinary or unexpected revelations. 
There are, however, some interesting points that should be dis- 
cussed: 
i. The general effect of the blades cn the meridional flow 
pattern is to increase the velocities near the shroud and 
decrease those near ihe hub. this may be considered as 
an improvement cf the meridional flow, for the most likely 
place tor separation is on the hub near station 13. Decreased 
velocities on the hub therefore lessen the likelihood of sepa- 
ration, 
2. the direction of the streamlines, given by 2, is changed 


very little by the action of the blades. 


AP 


ЕА at station 13 on the 
2 о 


3. The lowest static pressure, 
hub, is increased from -1.000 to -0.166 by the action of the 
blades. 


4. although relatively large changes of axial velocity occur 


within the blade region, the downstream eauilibrium values 


от 
differ very little (4.6 percent) from the axial velocity 
with no blades present. ilowever, tne tangential velocity 
far downstream is approximately proportional to the ra- 
dius, and is of the same orcer as the axial velocity: 
we = 1.146 at tip, = = 0.720 at root. Considerable 
diffusion would therefore be necessary to recover, as 
static pressure, ali of the kinetic energy added. 
5. Considerably more energy was added near the shroud 
than near the hub. ‘YWhis is unavoidable for radial blades, 
or, in fact, for any "practical" blades of tow aspect ratio. 
6. he total velocity increases discontinuously across 
the leading edge and it might Seem that th: pressure 
should decrease correspondingly. There is, however, a 
sufficient increase of total energy to cause the pressure 
to increase. his is because the tern: o rv is always 
greater than the term Lv in Zquation (124). 
7. The fact that the residuals are nearly constant on 
cylindrical surfaces extending downstream from the 
trailing edge of the blades means essentially that the 
tangential vorticity is nearly constant on these surfaces. 
rhis confirms, for this particular example at least, the 


linearizing assumption used by Marble (References 1 апа 
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4), when he assumed in the axial flow problem that the 
vorticity was transported unchanged on cylindrical sur- 
faces. It is surprising that this is alsc true in this 
region, bounded by cylindrical surfaces but located 


downstream of a complicated region of mixed flow. 





-%1. 
IM. AN ZSAMPLE А" SUBSsONIG. E LOW 
THROUGH AX ACTUATOR TOY 

an example of the motion of an incompressible fluid acted 
upon by an infinite system of blades was presented in Part VIII. ss 
regards rotationality effects, cormpressible fluid motion through a 
system of blades can be analyzed similerly, the main difference 
being the density ratio multiplying the vorticity term. Furthermore, 
it was shown in Part V-2 that tue flow in the region of tne blades, 
where tne vorticity is bound, is much less complicated than in the 
region downstream of the blades, when the vorticity is free. There- 
fore the present example, conceived in order to isolate the role - 
compressibility plays in altering the rotationality effects, is con- 
cerned only with the more complicated upstream and downstream 
regions. ihe blade region is therefore concentrated into an "ac- 
tuator disk' in which the tangential velocity jurnps discontinuously 
Írom zero to sorne finite value. It is assumed that the "blade sys- 
tem" is stationary, hence the total energy of the fluid is constant 
throughout the field. fhe matching conditions for the meridional 
flow at the discontinuity are the same as stated in Part V-D for 
leading edge discontinuities. 

The channel boundaries are concentric cylindrical surfaces, 
and at the entrance the axial velocity is constant and the radial and 


tangential velocities are zero. Lhe jump in tangential velocity is 
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equal in m agnivude to the inlet velocity cf the incompressible e: 





ample. 
Cor the downstream boundary condition the flow is independent of the 
«xial ccordinate. l'his problem is solved for two values of the inlet 
iach number (0 апа 0.2), the inlet mass-velocity and the tangential 
velocity being the sarae in both cases. 3 

Phe applicable equations, obtained from Zquations (63) and 
(65) with w= 0, (rv), =0, and хи = constant, are: 


Фф, = > 7 Fur, Pr (In e) y upstream 


der 


UA - 20] а (21)? downstream 
where rv(w) and E, ea are functions of W only and are evaluated 
at the downstream side of the discontinuity. 

“he boundary conditions are: Фф = constant on the cylindri- 
cal boundaries, -L/ = l on the upstream boundary, and Y, =0 
on the downstream boundary. For each approximation to rv( Y) 
the downstream boundary condition can be modified so that y is pre- 
Scribed, as in the previous example (Part VIII-E). 
A) Incompressible Flow Through the Actuator 
The difference equation corresponding to Zquation 125 for 


incompressible flow is: 


V, -(-21%, . + (ire )y, -4ф = О upstream ай 


1 (rv)? 2 downstream 
2 оқ 


S г?" 
ж-е 


The inlet veiocity is taken as 1000 and the boundary conditions are : 
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500 (9 - r’) on upstream boundary 
= Ооаг = 3 


4000 опг = 1 


й 


Ше € € € 
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0 оп downstream boundary 
rhe influence coefficients and the net points of the flow fieid are 
shown to scale in Figure 22. 

The first approximation to rv( y) was based cn irrotational 
flow within the same boundary. That is, it was assumed that the 
upstream value of y held throughout the field. The rv so obtained 
is shown in Figure 23. Before relaxation, the corresponding down- 
strearn boundary values were further improved by assuming half 
the change in Y cccurred unstream of the actuator, thus obtaining 
an approximate value of J at the actuator and consequently a new 
rv(Y). Tke boundary values corresponding to this last rv were 
then used for the first complete relaxation, after which rv( y) was 
found to be unchanged. The result of this complete relaxation was 
therefore the final solution of the problem. 

Values oí the axial, radial, and tangential velocities at sov- 
eral representative stations are shown in Figures 24, 25, and 26. 
Lines of constant static pressure are shown in Figure 27, Far down- 
stream of the Gctuator the axial velocity is much greater near the 
hub than near the shroud. The maximum value of the radial velocity 


occurs just be:ow the 'center' oí the channel. The pressure is 
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occur downstream of the actuator. 


Perhaps the most significant feature of the solution is that 


the radial velocity and the changes in axial velocity are approximately 





symmetrical about the actuator. Itis clear, however, that they are 
not exactly symmetrical since the forcing function is not symmetri- 
cal. One consequence of linearization (Reference 4) is that the mer- 
idional flow is symmetrical about the discontinuity. 

The solution cf this example will be discussed and compared 
with the compressible solution of the next section. 


в) Compressible flow Through the Actuator-Minjer = 0.2 





Che difference equation corresponding to Zquation (125) is: 
5 б 
d+ (IF Y eds (E lav 


Ly, (In г), / 5 á ` upstream (127)* 


=n loeh- e] SEU ве downstream 


The boundary conditions for the streamfunction* were essentially 


the same as for the previous example (IX-..). Cn the upstream 


boundary -2 Ls ew = 1000, and M = 0.2. Cther parameters on the 
2 


upstream bcundary were: 
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The tangential velocity jumped discontinuously to 1000 in the actuator. 
The influence coefficients and the net points of the field are 
the same as in the previous example (Figure 22). 

Phe first step in obtaining a solution was to estimate as ac- 
curately as possible, by any means whatever, the right side of the 
equation. This was done as follows: The vorticity term, ج‎ 
was very closely approximated by simply using the value from the 
incompressible solution just obtained. The extreme accuracy of 
this choice із indicated in Figure 28, where the initial and final val- 
ues of rv(y ) are plotted. The first estimation of the density term 
was likewise based on the irrotational velocities, but only values 
on the downstream boundary were calculated first. When this was 
done it was clear that the streamlines deflected less when compres- 
sibility was present, for the compressibility term, Й, (ind, 57. was 
negative and subtracted from the vorticity term, /-#($) 2 (го) 
which was positive. In addition f£] , which decreased according 
to the increase in total velocity downstream of the actuator, also 
acted to reduce the right side of the equation. However, the right 
е 
side of the equation was estimated on the downstream boundary using 
the density term just calculated, and the strearmfunction and the cor- 


responding density were calculated there. These values were treated 
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аз Comwnstresa. bcundary values for the first complete relaxation. 
in order to complete the estimetion of the residuals throughout the 
tie id the density changes calculated from the irrotational stream- 
lines were corrected by a factor chosen to give the correct density 
on the downstream boundary. in this way it was possible to make 
a rather accurate estimate of the right side of the equation, and 
the desired residuals fur the first complete relaxation. Drom the 
results of the first relaxation the densities and velocities, and 
finally new desired residuals, were obtained. fhe new desired 
residuals agreed very wel! with the previous ones except in regions 
where the Mach number was greatest. £mall changes of the stream- 
iunction brought the residuals into satisfactory agreement and led 
to the tinal solution, presented in Figures 28 through 32. 

the velocities are expressed in terms of the inlet axial 
velocity, rather than a velocity of sound, for comparison with the 
incompressible solution. Fhe Mach number based on the total ve- 
locity is shown in Figure 32. The velocity components are similar 


to those of the incompressibie solution, except oí course for the 


Giscontinuities at the actuator. 





С).. Discussion of tne Solutions 
in the two foregoing examples an otherwise axial flow is 
distorted by vorticity generated by a so-called actuator disk. [he 


Strength of the vorticity and the boundary conditions of the flow are 
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ae same in both examples, yet the solutions, although generally 
similar, ditfer distinctly in at least two ways: the total deflection 
of the flow brought about by the vorticity is decidedly less when 
the flow is compressible; the axial variations of the compressible 
fiow are concentrated more in the vicinity of the actuator, the 
concentration being rnore intense where the Mach number is higher. 

"he first point of difference was discussed in the previous 
section and is confirmed by Figure 33, where the streamlines of 
the two solutions are compared. A general conclusion set forth 
in Part VI is that the locai effect of vorticity, as regards deflec- 
tion of the streamsuríaces, reduces to zero as the "governing" 
“wach number increases to one, and then increases, but with орро- 
site sense, as the Mach number increases above one. The second 
point of difference could probably be predicted by the Prandtl- 
Glauert similarity transformation, applied te linearized equations. 
Its occurrence here is evident in Figure 32, and to some extent 
in Figures 25 and 30 where the maximum of the radial velocity 
for the compressible flow occurs nearer the center of the channel, 


indicating more rapid axial changes near the hub where the Mach 


number is larger. 
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This section provides examples of flow through actuator 
disks with maximum Mach numbers near one. Two examples, 
with the same inlet Mach number /0.555), but with actuator disks 
of different strength, are considered. In one example the maximum 
Mach number attained is one, in the other, 1.12. The actuator 
model is the same as described іа Part IX but the tangential ve- 
locity induced at the actuator is not constant there. Instead the 
tangential velocity is prescribed on the downstream boundary, and 
the jump at the actuator is not known until the problem is solved. 
With this prescription it is possible to replace the downstream 
boundary condition, = = 0, with the condition that the stream- 
function itself is prescribed, for the downstream conditions can 
be determined without knowing the intermediate flow. The limit- 
ing flow far downstream of the blades was discussed in Part V-7. 
For the particular case where the tangential velocity is constant, 
Figures 43 and 44 show the maximum Mach number as a function 
of the inlet Mach number and the tangential velocity. It was also 
Shown in Part V- that there is a maximum tangential velocity 
that can be imposed by an actuator disk and that this maximum 
corresponds to the choking condition, The limiting conditions 

for these examples are indicated by the two points shown in igure 
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|. fhe differentia} equation (Equation 125)) of the pre vious 


examples is applicable here but the definition of the streamfunc- 


tion* is slightly different. The ditference equation for both ex- 


amples is (from Zeuation (111) ): 
фун | фо т з (11 4, 
12 Los upstream (128) 
= 150, (2) ا‎ downstream 


3) Flow with Maximum Mach Number of C ne 
ihe following conditions are prescribed: 
inlet Mach number: MM = .555, м/ж ‚584 


Tangential velocity: у = 0 at inlet 


2. 
( = .¿l on downstream boundary 


and calculated frorn these are: 
Maximum axial velocity: em = .888 


М 
2 2? 
Maximum Mach number: = = 1 


The influence coefficients for each point in the net and the boundary 


values for Y are shown in Figure 22. 


The method oi solution was essentiaily the same as de- 


scribed in the previous example. In determining the approximate 


residuals for the first relaxation the value of the rotationality term 


/ rv |? 5 
2% =) obtained on the downstream boundary was assumed constant 


м” For these examples the streamfunction Y is defined as: 


2 e 
APA از‎ 4 
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ao С 


c» 
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on cylindrical surfaces extending upstream to the actuatcr. The 
density throughout the field was estimated using the known values 
of the density upstream and downstream and the approximate mage 
nitude of the density discontinuity at the actuator. The first de- 
sired residuals were obtained in this way and a complete relaxa- 
tion performed. The final solution then followed simple iteration 
in which the desired residuals were obtained from the solution of 
the previous step. After two complete relaxations the residuals 
were in good agreement except in the region where the Mach num- 
ber was largest. Minor modification then led to the final solution. 

he axiai and radial velocities are shown for representa- 
tive stations in Figures 35 and 36. The Mach number based on 
the total velocity is shown in Figure 37. 

his solution exhibits the properties already discussed іп 
the previous compressible example of “art IX, as noted in tha 
comparison o: the next section. 


В) Flow with 2 axiınum Wach Number of 1.12 











‚һе following conditions are prescribed: 
Inlet } ach number: M = „555, Wr = ‚584 


Tangential velocity v = O at inlet 


zZ 
(34 = .264 on downstrearr, boundary 


¥ 


and calculated from these are: 





MW. ,966 


Maximum axial velocity - 5 
+ 


Maxirnum axial Nach nurnber = 7 = .935 
Ж 
Maximum total Mach number = as = 1.12 

The influence coefficients for each point of the net and the boundary 
values for Y are shown in Figure 38. In addition the velocity, den- 
sity, etc., at the downstream boundary, are recorded in Table П. 

The residuals were first estimated using the meridicnal 
streamlines of the previous solution. A complete relaxation was 
performed and the iteration started. Unfortunately the densities 
obtained from the secend complete relaxation were imaginary and 
straightforward iteration could not be continued. To remedy this 
the axial mass-velocity, + Y, , was changed just enough to make 
the densities real. Then the iteration could be continued. It was 
observed that the numerical solution of this example was less ac- 
curate than the solution of the previous example, although the net 
points (Figures 34 and 38) were four times as dense. This is due 
to the nearly vertical tangent of the density curve (Figure 42) when 
the meridional velocity is nearly sonic, and to the fact that the ra- 
dial derivativa of the streamfunction could not be accurately deter- 
rnined on the inner radius where the large Mach numbers occur. 

The axiai and radial velocities are shown in figures 39 ana 
40, the Nach numbers in Figure 41. The solution exhibits the same 


properties as have already been described in the preceding examples. 





que 










it is surprising that the streamsurfaces are deflected less in this 
second example where the tangential velocity is greater. This can 
ое seen by comparing the downstreara values of the streamiunction 
in ;'igures 34 and 38, aud indicates that there is a certain tangen- 
tial deflection for which the defiection of the streamsuríaces is 


maximum. 











41. CONCIUSION 


. theory has been developed which permits the analysis of 
compressible flow in turbomachines having infinitely many blades 
when the governing velocity is subsonic. — examples, solved 
by the method oí finite differences, have been presented. The funda- 
mental idea underlying the theory is that the force field represent- 
ing the infinitely many blades is necessarily a '"pseudo-conservative' 
field. Because of this the three components of the field can be ex- 
pressed in terms of two functions, one describing the input of 
energy, the other the shape of the blades. The functions which 
must be prescribed and the boundary conditions which must be im- 
posed are then quite clear, and the heretofore more difficult direct 
problem becomes relatively easy. 

Two principal effects of compressibility were noted: the 
defiection of the streamsurfaces brought abcut by a given vorticity 
distribution decreased as the Mach numbers (subsonic) were in- 
creased; the streamwise variations of the flow became more con- 
centrated as the '.ach numbers were increased. 

Ihe foregoing theory is limited in that the fluid must be 
non-viscous and the number of blades must be infinite. In addition 
the application of the theory when the "governing" velocity is super- 
gonic is questionable. It is clear that non- viscous flow under the 


action of a finite number of blades must be completely understood 

















' vi$coud problem is attacked. Viscosity has no place in 
^ Ф . . з % LÀ . $ " e 
problems in which the number of blades is infinite, for the blades 
Y 


do not act as boundaries of the flow. 





The extension to me chines with a finite nuinber of blades 


could conceivably be carried out by three-dimensional finite differ- 


ence methods but the analysis would be very lengthy. A first order 


approximation to the pressure and total energy is indicated in the 


appendix. Per 


haps the least understood problems are those in which 


the governing velocity is transonic. The equations were developed 


and discussed generally in Part VI. + suggested approach might 


be to consider only the region in and upstream of the blades, where 


the vorticity effects are linear. Frescription of radial blades 


would allow further simplification by means of the transformation 


given in Zquation (121). Then a perturbation equation analogous to 


the transonic equations could be developed. Any analysis which 


will lead to a better understanding of the cushioning effect men- 


tioned in Part VI wouid be a worthy contribution. 
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^ first order approximation to the flow with a finite nuiabor 


of blades can be obtained frorn the solution with an infinite number 


rn 


сі blades by replacing the body force field acting in a sector between 


two blades with a pressure gradient force, thus essentiaily revers- 


ing the reasoning which first led to the assumption of infinite blades, 

Shen integration of the pressure gradient from one blade to another 

will give the circumferential pressure distribution between the bindes 
"he mean value of the pressure, as obtained from the anal- 


ysis with an infinite number of blades, is denoted by Ра“ Рр! 


denote the deviation of the total pressure from the mean: p- PtP! а 


¿he force field is then replaced - the pressure gradients as follows: 


/ 9, 
E / = г 2P | 
Ғ. = A кила -— o -r 
, 20" 


lf there are n blades distributed uniform ty around the axis, the 


second equation can be integrated from one blade to the next to give 


the nressure discontinuity across a blade: 


ar 
صد‎ = -}(S0) =A 

P 
Because lis independent of 9 , the tangential variation of the nres- 
sure is linear between blades. 


ihe deviation from the mean pressure 


can tuerefore be easily determined et all points in the region of the 
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vlades. The corresponding approximätion of the total energy, Y rv, 
cen be obtained from a form of Equation (31), using the velocities of 
the solution with an infinite number of blades. 

it is Significant that the equilibrium of the fluid has not been 
disturbed by the introduction of the pressure force for the force fie Иб, 
nor nas fulfillment of the continuity requirement нак altered. Ас- 
tually the only assuinption involved is that the tangential variation 
of the pressure is linear, i.e., A is independent of ©. This first 
order approximation therefore seems very reasonable, particularly 
when there are many blades, and may be satisfactory in other cases. 
Unfortunately exact solutions with a finite number of blades are not 


avcil-bie fcr determining the accuracy of the approximation. 
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ABLE I 


Shape Parameters for 
the Mixed Flow Example 


Shroud 
Radius 


3. 
За 
3. 
3. 


(м ім CQ) UJ WwW ly Go Ww 


LY Ww Ww w 


0900 
0000 
0000 
9600 
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.0072 
‚0227 
. 0489 


.0843 
.1250 
.1659 
‚2011 


.2273 
„2423 
. 2491 
. 2500 


. 2900 


Blade „паре 


Param.*, і, m 
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- .364 
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-.238 
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-.032 











фа ж pet ра iue ume ф- p 


Wn (lv [Do Ww 


.125 
. 25 
‚375 


„625 
.75 
„875 


-136- 


TABLE II 


Limiting Values Far Downstream for Transonic 
Actuator Zxample with Meee = 1.12 
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.9943 
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‚0626 
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ке ре рм (м 


ме 94 ме ме 


М ік: фр بيبط‎ фе 
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„0691 


‚0869 
. 1033 
.1198 
‚1351 
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.1774 
.1904 
.2029 


tor 
Y^ ап 
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.21456 
‚19310 
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‚16054 
‚14793 
‚13709 


‚12767 
. 11942 
‚11214 
. 10566 


‚09985 
‚09464 
. 06993 
.08564 
.08173 


. 9658 
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„9027 
‚8744 
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1132 .1290 .1458 .1655 .1882 .2022 .2173 .2350 .2550 
‚1764 .1936 .2116 .2304 .2500 .2601 .2704 .2809 .2916 
‚1605 .1782 .1967 .2160 .2363 .2465 .2574 .2654 .2796 
,1445 .1627 .1817 .2017 .2226 .2334 .2444 .2559 .2675 
‚1286 .1473 .1668 .1845 .2089 .2201 .2314 .2434 . 2550 
‚1126 .1318 .1518 .1729 .1952 .2067 .2185 .2309 .2436 
„0967 .1164 .1369 .1586 .1814 .1934 .2055 .2184 .2315 
,0807 .1009 .1220 .1437 .1677 .1800 .1925 .2059 .2195 
‚0648 .0855 .1070 .1298 .1540 .1667 .1795 .1934 .2075 
„0488 .0700 .0921 .1155 .1403 .1533 .1665 .1809 .1955 
.0329 .0546 .0771 .1011 .1266 .1400 .1535 .1685 .1835 
.0169 .0391 .0622 .0867 .1129 .1266 .1405 .1560 .1715 
¿0010 .0237 .0473 .0724 .0992 .1133 .1276 .1435 .1595 

,0092 .0323 .0590 .0855 .0999 .1146 .1310 .1475 
.0174 .0436 .0717 .0866 .1036 .1185 .1354 
.0024 .0293 .0580 .0733 .0886 .1060 .1234 

.0149  .0443 .0599 .0756 .0935 .1114 
.0005 .0306 .0406 .0626 .0810 .0994 
.0169 .0332 .0466 .0685 .0874 
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TABLE Ill (Cont'd) 
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„20 
.24 
.26 
. 20 


.30 
soe 
. 34 
. 36 
. 38 


‚40 
742 
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‚46 
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.2800 .3090  .3300 .3525 .3836 4550  .5350 .5809 
3025. .3136 — .3192 .3249 .3306 .3345 .3306 .3249 
2911 .3028 .3089 .3150 .3213 .3264 .3237 .3186 
2796 .2920 .2985 .3051 .3120 .3184 .3169 .3124 
‚2682 .2813 .2882 .2952 .3027 .3103 .3100 .3061 
‚2568 .2705 .2778 .2854 .2934 .3023 .3032 .2998 
-2453 .2597 .2675 .2755 .2841 .2942 .2963 .2936 
-2339 .2489 .2572 .2656 .2748 .2862 .2895 .2873 
‚2225 .2382 .2468 .2557 .2655 .2781 .2826 .2811 
2111 .2274 — .2365 .2458 .2562 .2701 .2758 .2748 
1996 .2166 | .2261 .2359 .2969 .2620 .2689 .2685 
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1425 .1627  .1744 .1865 .2004 .2218 .2347 .2372 
1310 .1522 .1641 .1766 .1911 .2137 .2278 .2309 
1195 .1412 .1538 .1668 .1818 .2057 .2210 .2247 
1082 .1304 .1434 .1569 .1725 .1976 .2141 .2184 
0968 .1196 .1331 .1470 .1632 .1896 .2073 22121 
0853 .1089 .1227 .1371 .1539 .1815 .2004 .2059 
0739 .0981  .1124 .1272 .1446 .1735 .1936 .1996 
0625 .0873  .1021 .1173 .1353 .1654 .1867 .1934 
МИЫ Ой Cor Ме 150 04 11909. „0671 
0396 .0658 .0814 .0976 .1167 .1493 .1730 .1808 
0282 .0550 .0710 .0877 .1074 .1413 .1662 .1746 
‚0167 .0440 .0607 .0778 .0981 .1332 .1593 .1683 
‚0053 .0334 .0504 .0679 .0889 .1251 .1525 .1620 
0222. .0400 .0580 .0796 .1171 .1456 .1558 

0119. .0297 .0481 .0703 .1090 .1398 .1495 

001: (0193 .0383 .0610 .1010 .1319 .1432 

0090 .0284 .0517 .0929 .1251 .1370 
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шо چ‎ da UI СУ 


DI Po 














Inge ‚6450 .6955 9802 9074 
‚3136 .3625 2916 .2809 .2704 .2601 .250C .2304 
.3081 .2975 2871 .2757 .2666  .2565 .2467 .2276 
‚3026 .2925 2825 .2726 .3M27 .2530 .2434 .2247 
.2971 .2875 2780 .2684 .2569 .2494 .2402 .¿219 
‚2916 .2826 2735 .2643 .2550 .2459 .2369 .219%6 
.2960 .2776 2639 .2601 .2512 .2423 .2336 .2162 
‚2805 .2726 2644 .2560 .2474 .2355 .2301 .2133 
‚2750 .2676 2599 .2518 .2435 .2352 .2270 .2165 
‚2695 .2626 2553 .2477 .2397 .2317 .2238 .2077 
2640 .2576 2508 .2435 .2358 .2281 .2205 .2048 
„25865 .2527 .2463 .2393 .2320 246 .2172 .2020 
„2920 .2477 .:417 .2552 .2282 210 .2139 .1991 
‚2475 .2427 2372 .2310 .224% ,2175 .2106 .1963 
‚2419 „2377 2327 .2208 .2205 .2139 .2074  .1935 
‚2364 .2327 2282 .2227 .2166 .2103 .2041 .1°06 
2309 .2277 222 2185 .2128 .2068 .2008 .1878 
.2254 .2228 2101 .2144 .2090 .2032 .1975 .1845 
,2199. .2178 .2146 .2102 .2051 .1997 .1942  .1621 
.2144 : ,2128 2100 .2061 .2013 .1961 .1910 .1792 
„2089 .2078 2055 .2010 .1974 .1926 .1877 .1764 
.2034 .2028  .2010 .1977 .1936 .1890 .1844 .1736 
.1978  .1978 .1967 .1936 .1896 .1555 .1811 .1707 
‚1923 .192 .1919 .1894 .1859 .1819 .1778 .1679 
‚1868  .1679 .1874 .1353 .1821 .1784 .1746 .1658 
.1813 .1822 .1828 .1811 .1782 .1782 .1713 .1622 
„1758 979 .1783 .176% .1743 .1712 .1673 .1593 
‚1703 .1725 |  .1732 .1728 .1706 .1677 .1643 .1565 
‚1648 .1979 .1692 .1686 .1667 .1641 .1614 .1531 
‚1593 .1629 .1647  .1645 .1629 .1606 .1582  .1508 
.1538 .1580 .1602 .i633 .1590 .1570 .1549  .1480 
.l482 .1530 .1556 .156. .1532  .1535 .1316 .1451 






































